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SYNOPSIS 



Supersymmetry remains the best motivated extension of the Standard Model, since it 
stabilizes the gauge hierarchy against quadratically divergent radiative corrections. How- 
ever, this argument by itself does not yet tell us very much about the way supersymmetry 
is realized in Nature. In particular, it tells us nothing about the mechanism that breaks 
supersymmetry, except that sparticles masses should not greatly exceed 1 TeV (with the 
possible exception of first and second generation sfermions; this possibility will not be pur- 
sued here). One can be completely "agnostic" and simply parameterize SUSY breaking by 
soft terms. This approach has largely been taken by the first and second subgroup of this 
working group. Alternatively, one can explore the predictions of specific models, as has 
been done by subgroups three (on gauge mediated SUSY breaking) and four (on anomaly 
mediated SUSY breaking). 

An additional complication arises because supersymmetric models allow new, gauge- 
invariant interactions that have no equivalent in the standard model; this is true even if one 
sticks to the minimal particle content. These interactions break R— parity. Some implications 
of this have been investigated in subgroup two. 

Some aspects of SUSY phenomenology have already reached quite a high level of maturity. 
Radiative corrections have been computed for total cross sections and partial widths; their 
impact on kinematic distributions is now being explored. Finite width effects formally also 
belong to higher orders in perturbation theory. Finally, in some cases it is important to 
keep spin correlations through the entire sparticle production and decay process. Some 
investigations along these lines are described in contributions from subgroup one; they are 
(possibly with minor modifications) relevant for all SUSY models. 

Another area of current interest is the determination of soft breaking parameters, in 
particular of sparticle masses. Here a realistic assessment of the capabilities of various 
experiments requires detailed Monte Carlo simulations. Several studies of this kind are 
described in contributions from subgroups one, two and four. 

Manuel Drees (for the Working Group). 
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SUSY particle production at hadron colliders 

M. SPIRAQ 



Abstract 

The determination of the full SUSY QCD corrections to the production of squarks, 
gluinos and gauginos at hadron colliders is reviewed. The NLO corrections stabilize the 
theoretical predictions of the various production cross sections significantly and lead to 
sizeable enhancements of the most relevant cross sections. We discuss the phenomeno- 
logical consequences of the results on present and future experimental analyses. 

1 Introduction 

The search for Higgs bosons and supersymmetric particles is among the most important en- 
deavors of present and future high energy physics. The novel colored particles, squarks and 
gluinos, and the weakly interacting gauginos can be searched for at the upgraded Tevatron, 
a pp collider with a cm. energy of 2 TeV, and the LHC, a pp collider with a cm. energy 
of 14 TeV. Until now the search at the Tevatron has set the most stringent bounds on the 
colored SUSY particle masses. At the 95% CL, gluinos have to be heavier than about 180 
GeV, while squarks with masses below about 180 GeV have been excluded for gluino masses 
below ~ 300 GeV [0. Stops, the scalar superpartners of the top quark, have been excluded 
in a significant part of the MSSM parameter space with mass less than about 80 GeV by 
the LEP and Tevatron experiments Q. Finally charginos with masses below about 90 GeV 
have been excluded by the LEP experiments, while the present search at the Tevatron is 
sensitive to chargino masses of about 60-80 GeV with a strong dependence on the specific 
model PJ. Due to the negative search at LEP2 the lightest neutralino Xi nas to be heavier 
than about 30 GeV in the context of SUGRA models JIJ. In the i?-parity-conserving MSSM, 
supersymmetric particles can only be produced in pairs. All supersymmetric particles will 
decay to the lightest supersymmetric particle (LSP), which is most likely to be a neutralino, 
stable thanks to conserved i?-parity. Thus the final signatures for the production of su- 
persymmetric particles will mainly be jets, charged leptons and missing transverse energy, 
which is carried away by neutrinos and the invisible neutral LSP. 

In Section 2 we shall summarize the details of the calculation of the NLO QCD correc- 
tions, as described in Refs. |2|-§J for the case of qq production. The evaluation of the full 
SUSY QCD corrections splits into two pieces, the virtual corrections, generated by virtual 
particle exchanges, and the real corrections, which originate from gluon radiation and the 
corresponding crossed processes with three-particle final states. 

In Section 3 we shall consider the production of squarks and gluinos except stops 0. We 
assume the light-flavored squarks to be mass degenerate, which is a reasonable approximation 
for all squark flavors except stops, while the light quarks (u, d, s, c, b) will be treated as 

*Heisenberg Fellow 



6 



massless particles. The production of stop pairs requires the inclusion of mass splitting and 
mixing effects || and will be investigated in Section 4. In Section 5 we will summarize 
the results for the production of charginos and neutralinos at NLO Qj. The calculation 
of the LO cross sections has been performed a long time ago [||. Since in most of the 
cases the [unphysical] scale dependence of the LO quantities amounts up to about 50%, the 
determination of the NLO corrections is necessary in order to gain a reliable theoretical 
prediction, which can be used in present and future experimental analyses. 

2 SUSY QCD corrections 
2.1 Virtual corrections 



g nrsmywy 
i 




g ^SLQSLSJU 



Figure 1: Typical diagrams of the virtual corrections. 



The one-loop virtual corrections are built up by gluon, gluino, quark and squark exchange 
contributions [see Fig. They have to be contracted with the LO matrix elements. The 
calculation of the one-loop contributions has been performed in dimensional regularization, 
leading to the extraction of ultraviolet, infrared and collinear singularities as poles in e = 
(4 — n)/2. For the chiral 75 coupling we have used the naive scheme, which is well justified 
in the present analysis at the one-loop levelQ. We have explicitly checked that after summing 
all virtual corrections no quadratic divergences are left over, in accordance with the general 
property of supersymmetric theories. The renormalization has been performed by identifying 
the squark and gluino masses with their pole masses, and defining the strong coupling in the 
MS scheme including five light flavors in the corresponding (3 function. The massive particles, 
i.e. squarks, gluinos and top quarks, have been decoupled by subtracting their contribution 
at vanishing momentum transfer ||. In dimensional regularization, there is a mismatch 
between the gluonic degrees of freedom [d.o.f. = n — 2] and those of the gluino [d.o.f. = 2], 
so that SUSY is explicitly broken. In order to restore SUSY in the physical observables in 
the massless limit, an additional finite counter-term is required for the renormalization of 
the novel qgq vertex 0. 

T We have explicitly checked that the results obtained with a consistent 75 scheme are identical to the one 
with the naive scheme. 
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2.2 Real corrections 




9 



9 



9 



Figure 2: Typical diagrams of the real corrections. 



The real corrections originate from the radiation of a gluon in all possible ways and 
from the crossed processes by interchanging the gluon of the final state against a light 
(anti)quark in the initial state. The phase-space integration of the final-state particles has 
been performed in n = 4 — 2e dimensions, leading to the extraction of infrared and collinear 
singularities as poles in e. After evaluating all angular integrals and adding the virtual 
and real corrections, the infrared singularities cancel. The left-over collinear singularities 
are universal and are absorbed in the renormalization of the parton densities at NLO. We 
defined the parton densities in the conventional MS scheme including five light flavors, i.e. 
the squark, gluino and top quark contributions are not included in the mass factorization. 
Finally we end up with an ultraviolet, infrared and collinear finite partonic cross section. 

However, there is an additional class of physical singularities, which have to be regu- 
larized. In the second diagram of Fig. ^| an intermediate qg* state is produced, before the 
[off-shell] gluino splits into a qq pair. If the gluino mass is larger than the common squark 
mass, and the partonic cm. energy is larger than the sum of the squark and gluino masses, 
the intermediate gluino can be produced on its mass-shell. Thus the real corrections to qq 
production contain a contribution of qg production. The residue of this part corresponds to 
qg production with the subsequent gluino decay g — > qq, which is already contained in the 
LO cross section of qg pair production, including all final-state cascade decays. This term 
has to be subtracted in order to derive a well-defined production cross section. Analogous 
subtractions emerge in all reactions: if the gluino mass is larger than the squark mass, the 
contributions from g — > qq, qq have to be subtracted, and in the reverse case the contributions 
of squark decays into gluinos have to subtracted. 



Squarks and gluinos can be produced via pp,pp — > QQ,QQ,Qg,gg a t hadron colliders. The 
hadronic squark and gluino production cross sections can be obtained from the partonic 
ones by convolution with the corresponding parton densities. We have performed the nu- 
merical analysis for the upgraded Tevatron and the LHC. For the natural renormaliza- 
tion/factorization scale choice Q = m, where m denotes the average mass of the final-state 
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SUSY particles, the SUSY QCD corrections are large and positive, increasing the total cross 
sections by 10-90% 0. This is shown in Fig. || where the K factors, defined as the ratios 
of the NLO and LO cross sections, are presented as a function of the corresponding SUSY 
particle mass for the LHC. We have investigated the residual scale dependence in LO and 
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Figure 3: K factors of the different squark and gluino production cross sections at the LHC 
[yfs = 14 TeV]. Parton densities: GRV94 with Q = m. Top mass: m t = 175 GeV. 

NLO, which is presented in Fig. |j. The inclusion of the NLO corrections reduces the LO 
scale dependence by a factor 3-4 and reaches a typical level of ~ 15%, which serves as an 
estimate of the remaining theoretical uncertainty. Moreover, the dependence on different 
sets of parton densities is rather weak and leads to an additional uncertainty of ~ 15%. In 
order to quantify the effect of the NLO corrections on the search for squarks and gluinos 
at hadron colliders, we have extracted the SUSY particle masses corresponding to several 
fixed values of the production cross sections. These masses are increased by 10-30 GeV at 
the Tevatron and 10-50 GeV at the LHC, thus enhancing the present and future bounds on 
the squark and gluino masses significantly. Finally we have evaluated the QCD-corrected 
transverse-momentum and rapidity distributions for all different processes. As can be in- 
ferred from Fig. |5|, the modification of the normalized distributions in NLO compared to 
LO is less than about 15% for the transverse- momentum distributions and much less for the 
rapidity distributions. Thus it is a sufficient approximation to rescale the LO distributions 
uniformly by the K factors of the total cross sections. 
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Figure 4: Sca/e dependence of the total squark and gluino production cross sections at the 
LHC in LO and NLO. Parton densities: GRV94 (solid), CTEQ3 (dashed) and MRS(A') 
(dotted); mass parameters: rriq = 600 GeV, m g = 500 GeV and m t = 175 GeV. 



4 Stop Pair Production 

At LO only pairs of t\ or pairs of I2 can be produced at hadron colliders. Mixed t\t2 pair 
production is only possible at NLO and beyond. However, we have estimated that mixed 
stop pair production is completely suppressed by several orders of magnitude and can thus 
safely be neglected The evaluation of the QCD corrections proceeds along the same lines 
as in the case of squarks and gluinos. The strong coupling and the parton densities have been 
defined in the MS scheme with 5 light flavors contributing to their scale dependences, while 
the stop masses are renormalized on-shell. The QCD corrections increase the cross sections 
by up to about 40% [see Fig. |6| and thus lead to an increase of the extracted stop masses 
from the measurement of the total cross section. Moreover, as in the squark/gluino case the 
scale dependence is strongly reduced [see Fig. |7j and yields an estimate of about 15% of the 
remaining theoretical uncertainty at NLO. At NLO the virtual corrections depend on the 
stop mixing angle, the squark, gluino and second stop masses. However, it turns out that 
these dependences are very weak and can safely be neglected as can be inferred from Fig. [8[ 
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Figure 5: Normalized transverse-momentum and rapidity distributions of pp — > qq + X at 
the upgraded Tevatron [yfs = 2 TeV] in LO (dotted) and NLO (solid). Parton densities: 
CTEQ4L (LO) and CTEQ4M (NLO) with Q = m; mass parameters: nig = 250 GeV, 

175 GeV. 
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Figure 6: K factor of the light stop production cross sections at the upgraded Tevatron 
[y/s = 2 TeV]. Parton densities: CTEQ4L (LO) and CTEQ4M (NLO) with Q = m h . Top 
mass: m t = 175 GeV. 



5 Chargino and Neutralino Production 

The production cross sections of charginos and neutralinos depend on several MSSM param- 
eters, i.e. Mi,M2,/x and tan/3 at LO 0. The cross sections are sizeable for chargino/neu- 
tralino masses below about 100 GeV at the upgraded Tevatron and less than about 200 
GeV at the LHC. Due to the strong dependence on the MSSM parameters the extracted 
bounds on the chargino and neutralino masses depend on the specific region in the MSSM 
parameter space ffl. The outline of the determination of the QCD corrections is analogous 
to the previous cases of squarks, gluinos and stops. The resonance contributions due to 
gq — > XiQ with q q\j have to be subtracted in order to avoid double counting with the 
associated production of gauginos and strongly interacting squarks and gluinos. The parton 
densities have been defined with 5 light flavors contributing to their scale evolution in the 
MS scheme, while the t-channel squark masses have been renormalized on-shell. The QCD 
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Figure 7: Scale dependence of the total light stop production cross sections at the Tevatron 
and LHC in LO and NLO. Parton densities: CTEQ4L (LO) and CTEQ4M (NLO); Top 
mass: m t = 175 GeV. 

corrections enhance the production cross sections of charginos and neutralinos by about 10- 
40% [see Fig. |J. The LO scale dependence is reduced to about 10% at NLO [see Fig. |], 
which signalizes a significant stabilization of the thoretical prediction for the production 
cross sections The dependence of the chargino/neutralino production cross sections on 
the specific set of parton densities ranges at about 15%. 

6 Conclusions 

In this work we have reviewed the status of SUSY particle production at hadron colliders at 
NLO. Most QCD corrections to the production processes are known, thus yielding a nearly 
complete theoretical status. There are especially large QCD corrections to the production 
of gluinos, which significantly increase the extracted bounds on the gluino mass from the 
negative search for these particles at the Tevatron. In all processes, which are known at 
NLO, the theoretical uncertainties are reduced to about 15%, which should be sufficient for 
the upgraded Tevatron and the LHC|. 
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[dashed] and NLO [solid]. The thickness of the NLO curves represents the dependence of the 
cross sections on the stop mixing angle and the squark and gluino masses. The shaded NLO 
bands indicate the theoretical uncertainties due to the scale dependence within mj/ '2 < \i < 
2m- t . Parton densities: CTEQ4L (LO) and CTEQ4M (NLO); Top mass: m t = 175 GeV. 
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Comparison of exact matrix element calculations with ISAJET 
and PYTHIA in case of degenerate spectrum in the MSSM 

S. ABDULLIN, V. ILYIN and T. KON 
Abstract 

We calculate the transverse momentum distribution of additional partons from 
VP * 99 + X with nig = 1.2 TeV at y / i=14 TeV. The calculations are carried out 
with exact 2^3 matrix element generators CompHEP and GRACE, and widely used 
physics generators ISAJET and PYTHIA. In the latter cases only 2^2 subprocesses 
are included and additional partons originate mainly from initial and final states radi- 
ation. The results obtained by CompHEP and GRACE are in perfect agreement in the 
whole region of px (10 - 2000 GeV) of interest. The comparison with p? distributions 
produced by ISAJET and PYTHIA does not show significant difference. ISAJET tends 
to overshoot for px values in the region of 200-900 GeV, while PYTHIA's distribution 
has quite a good agreement with exact ME calculations. 

One of the main purposes of the LHC collider is to search for physics beyond the Stan- 
dard Model (SM), in particular to look for superpartners of ordinary particles expected in 
SUperSYmmetric extensions of the SM (SUSY). SUSY, if it exists, is expected to reveal 
itself at the LHC firstly via an excess of (multilepton+) multijet + E™ lss final states com- 
pared to Standard Model (SM) expectations The bulk of these studies was carried 
out in the framework of the Minimal SuperSymmetric Model (MSSM), or the Supergravity- 
inspired subset of MSSM, minimal SUGRA (mSUGRA) with universal gaugino mass and 
corresponding mass hierarchy. 

Unfortunately even in the framework of the MSSM one can expect scenarios with non- 
universal gaugino mass leading to some "pathological" SUSY cases which are difficult to 
observe due to a "degenerate" mass spectrum [[|. For instance, it is quite obvious that there 
are no visible objects in the SUSY event if m^o « m g m nig = Msusy- In this case one can 
rely mainly on observation of jet activity which is additional to sparticle pair production. 
The idea is that additional high p? recoil partons can provide a clear E™ lss +jet(s) signature 
to observe the excess of SUSY events over SM expectations. 

Our preliminary estimates of the MSSM signal generated with ISAJET [|J] show that the 
limit of visibility (at the 5 a level) in the MSSM is about 1.2 - 1.3 TeV with integrated 
luminosity of 100 fb -1 for m^o « m g <C rriq, thus much smaller than the squark-gluino mass 
reach (> 2.5 TeV) in mSUGRA with the same integrated luminosity ||. 

It is well known that data bases in ISAJET and PYTHIA ||] contain matrix elements 
only for 2 — > 2 subprocesses. So, additional partons can be produced only via initial and final 
state radiation (ISR/FSR parton showering). If one discusses production of hard partons, 
this mechanism is an approximation with a priori unknown precision. From a general point 
of view this approximation should underestimate the result obtained within the perturbative 
QCD theory through evaluation of the complete set of relevant Feynman diagrams. Thus, 
one might assume that there will be a substantial number of events with high pt of the heavy 
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SUSY system recoiling against observable jet(s), larger than predicted by the parton shower- 
ing mechanism in ISAJET and PYTHIA. We use the packages CompHEP @ and GRACE 
to evaluate 2 — > 3 matrix elements in leading order and calculate the distributions. Finally 
we compare these results with ones obtained with ISAJET and PYTHIA. 

Calculational details 

To make a coherent calculation we choose a set of parameters as follows: 

• gluino mass m g = 1.2 TeV; 

• common squark mass m g = 100 TeV (10 TeV for PYTHIA) to neglect squark contri- 
butions; 

• CTEQ4M set of structure functions; 

• Q 2 = Am 2 - ( in PYTHIA Q 2 = (m 2 ±1 + m 2 ±2 )/2); 

• NLO running strong coupling with normalization ots{Mz) = 0.118, so as(2400 GeV) = 
0.0817. 

We use ISAJET 7.44 and PYTHIA 6.125 with multiple interactions switched off. We treat 
the jet in ISAJET/PYTHIA not as the single hardest final state parton, but use the final 
state parton list (excluding the two gluinos) as input to an iterative cone jet finder algorithm 
(Rcone = 0.5) to find a group of partons forming parton- level jets. 

In CompHEP and GRACE we calculate exact tree level matrix elements for the subpro- 
cesses : 

gg^gg + g; gq(Q) -> gg + q(Q); qQ^gg + g; 

which were convoluted with the CTEQ4m parton distributions to get LHC cross sections 
and distributions. The q (Q) terms stand for quarks (antiquarks) of the first five flavours: 
u, d, c, s and b. 

We ask the final (additional) parton, gluon or quark, to radiate within the typical ATLAS 
and CMS rapidity, \r)j\ < 4.5, and to be hard enough, E 3 T > 10 GeV. No cuts are applied on 
the gluino pair. 

Results and Conclusions 

The cross section of the process pp — > gg + jet within cuts on the jet turns out to be rather 
high, about 120 fb, with the following contributions from different channels : 

a gg = 78.9 fb; a gq{Q) = 21.9 fb; a qQ = 20.0 fb. 

We found that the shapes of the pr distributions are quite different in these three cases, 
see Fig.|l|. Indeed, one can expect that the qQ channel produces less events than the gg 
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contribution, due to the dominance of the gluon component in the proton for small x. 
However, one can see that the gq{Q) curve is higher than the gluon-gluon one for > 400 
GeV, being much smaller even than qQ at the soft end. Note that the set of topologies of 
Feynman diagrams is the same in all three channels. So, the observed difference connects 
surely with different interplay of contributions of vector and spinor virtual particles. 

In Fig.|2] the p? distribution for pp —> gg + jet is given as a sum of the three channels 
discussed above. The CompHEP/GRACE result is represented by solid histograms, IS A JET 
by dashed histograms and PYTHIA by dot-dashed ones. One can see that the PYTHIA 
distribution underestimates the exact ME result smoothly. Of course this (almost constant) 
difference can be explained (corrected) by a different setting of some QCD switches. At the 
same time one can see a more pronounced deviation from the exact ME in the shape of the 
ISAJET curve. In particular at moderate p? (~ 200 — 900 GeV) ISAJET overestimates the 
exact ME result. One can also see that for very high transverse momenta (> 1200 GeV for 
ISAJET and 1400 GeV for PYTHIA) the parton showering approximation is too crude. 

Coming back to the initial idea to look at additional jet activity with high p? recoiling 
against the heavy SUSY system one can conclude that parton shower approximation works 
well enough up to very high values of transverse momenta. Therefore, unfortunately one can 
not find room for sizable improvements of parton shower predictions by using exact matrix 
elements in the process discussed. 

Some interesting discrepancy is observed in the region of relatively small transverse mo- 
menta, pt < 150 GeV, where ISAJET and PYTHIA produce much weaker (by a factor of 
2-3) jet activity than CompHEP/GRACE. However, this region is sensitive to higher order 
corrections (in particular due to resummation of large logarithms of p^jm-^), thus a more 
careful analysis is necessary to get reliable conclusions. 
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Figure 1: Contribution of different channels to the pt distribution of the additional jet in 
PP 99 + X, calculated with CompHEP and GRACE. 
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Figure 2: distribution of the additional jet in pp — > gg + X, calculated with different 
packages. 
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The finite width effect on neutralino production 

T. KON, Y. KURIHARA, M. KURODA and K. ODAGIRI 

It has been proposed in || that neutralino mass differences can be determined fairly 
accurately through the /x-pair invariant mass distribution da/dMyy of the process 

PP -> xtxl x -> xt(pfi)X -> XiViXiriXi C 1 ) 

as the distribution exhibits a sharp decrease at its kinematic endpoint. When, however, the 
width of the neutralino is taken into account, the smearing of the cross section might make 
this observation obsolete and less appealing for the determination of the neutralino mass. 
Therefore, it is important to know to what extent the width of neutralino makes the mass 
determination less accurate. 

We have investigated this problem by comparing the muon-pair invariant mass distribu- 
tion in the following three cases. 

(case 1) Zero width approximation for Only those diagrams which proceed via the chain 
decay (1) are taken into account. 

(case 2) Same as (1) but finite width of x® * s used. Only those diagrams that proceed via 
resonances as shown in (1) are considered. 

(case 3) Finite width of x° as we U as the entire set of diagrams (165 diagrams in unitary 
gauge) that creates the final state XiXi^V are taken into account. The diagrams which are 
considered in case 3 but not in case 2 constitute the background to this process. 

In the numerical evaluation of the muon pair distribution, we have used the following set 
of SUSY parameters, 

tan/? = 12, M 2 = 200GeV, /i = -500GW, (2) 

which results in the following masses for charginos and neutralinos, 

m _o = (97,197,507,511), (3) 
m-+ = (197,514). (4) 

In particular, in the present set of parameters we are considering, Xi is almost purely bino 
(the SUSY partner of the U(l) gauge boson, By), while x° i s almost purely neutral wino 
(the SUSY partner of the neutral SU{2) L gauge-boson, W%). As for the slepton masses, we 
have used the common mass m| =150 GeV. 

Assuming that other SUSY particles are heavy so that x° decays dominantly into ££, we 
have evaluated the width of x® as 

T^o = 1.25GeV. (5) 
The process is described by the subprocess 

QQ xtxl -> Xi(vfi) -> XiKx.iV>)- (6) 
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For the quark distribution functions, we have used QTEC-3D [0. 

The numerical evaluation was performed by HERWIG H in case 1 and case 2, while the 
case 3 was calculated by GRACE [|]. 

Results are shown in Fig. 1(a) and (b), where the distributions of the invariant mass 
in the second neutralino decay are shown. In Fig. 1(a), the dashed line represents the 
zero width approximation (case 1), the dotted line represents the distribution for the finite 
X® width but with only the resonance diagrams (case 2), while the solid line represents the 
distribution for the finite x° width and full set of diagrams which give the final states xfxi^f 1 
(case 3). We show a detailed comparison between case 2 and case 3 around the end point 
region of M MAt in Fig. 1(b). As one sees from the figures, the edge of the distribution at the 
end point of the phase space remains, although it is broadened by the finite width effect. As 
a result, the exact position of the end point acquires an ambiguity of about 5 GeV, although 
the precision of mass determination, which may additionally make use of the shape of the 
distribution, can presumably be made much less than this. A full analysis will need to take 
into account the experimental resolutions and signal selection efficiency, and this is beyond 
the scope of our study. The total cross section of the signal process (1) is 0.015 pb. 

Acknowledgements One of the authors, K.O., would like to thank Drs. J. Kanzaki 
and P. Richardson for their computing help. 
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Figure 1: Distribution of the invariant mass in the second neutralino decay from the 
X1X2 associated production. In (a) a comparison is made among the full calculation using 
GRACE (continuous line, case 3), the calculation with HERWIG assuming a narrow width and 
an isotropic decay (dashed line, case 1), and another calculation using HERWIG incorporating 
a finite X2 width (dotted line, case 2). (b) corresponds to the detailed comparison between 
case 2 and case 3 around the end point region of M w . 
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Width effects in slepton production e + e — > j^^j^p 

H.-U. MARTYN 

Abstract 

A case study will be presented to investigate width effects in the precise determi- 
nation of slepton masses at the e + e~ Tesla Linear Collider. 

1 Introduction 

If supersymmetry will be discovered in nature a precise measurement of the particle spectrum 
will be very important in order to determine the underlying theory. The potential of the 
proposed Tesla Linear Collider [T[| with its high luminosity and polarization of both e ± 
beams will allow to obtain particle masses with an accuracy of 10~ 3 or better [Q]. At such a 
precision width effects of primary and secondary particles may become non- negligible. 

The present case study is based on a particular i?-parity conserving mSUGRA scenario, 
also investigated in the Ecfa/Desy Study |§, with parameters m = 100 GeV, m 1( / 2 = 
200 GeV, A = GeV, tan/5 = 3 and sgn(/i) > 0. The particle spectrum is shown in 
fig. [p. Typical decay widths of the scalar leptons are expected to be T ~ 0.3 — 0.5 GeV, 
while the widths of the light gauginos, decaying into 3-body final states, are (experimentally) 
negligible. 




XT 



Figure 1: Mass spectrum and decay modes of sleptons and light gauginos 
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This note presents, as an example, a simulation of right scalar muon production 




(1) 



The analysis is based on the methods and techniques described in a comprehensive study of 
the same Susy spectrum The detector concept, acceptances and resolutions are taken 
from the Tesla Conceptual Design Report |L|. Events are generated with the Monte Carlo 
program Pythia 6.115 ||, which includes the width of supersymmetric particles as well 
as QED radiation and beamstrahlung ||. It is assumed that both beams are polarized, 
right-handed electrons to a degree of V e - = 0.80 and left-handed positrons by V e + = 0.60. 
A proper choice of polarizations increases the cross section by a factor of ~ 3 ancl reduces 
the background substantially, e.g. by more than an order of magnitude for Standard Model 
processes. 

2 Mass determinations 

Scalar muons fiR are produced in pairs via s channel 7 and Z exchange and decay into an 
ordinary muon and a stable neutralino x? (LSP), which escapes detection. The experimental 
signatures are two acoplanar muons in the final state with large missing energy and nothing 
else in the detector. Simple selection criteria [|IJ (essentially cuts on acollinearity angle and 
missing energy) suppress background from W + W~ pairs and cascade decays of higher mass 
Susy particles and result in detection efficiencies around ~ 70%. 

Two methods to determine the mass of /2r will be discussed: (i) a threshold scan of the 
pair production cross section, and (ii) a measurement of the energy spectrum of the decay 
muons, which simultaneously constrains the mass of the primary smuon and the secondary 
neutralino. The particle mass parameters given by the chosen Susy model are mp -R = 
132.0 GeV, T^ R = 0.310 GeV and m x? = 71.9 GeV. 

2.1 Threshold scan 

Cross section measurements close to production threshold are relatively simple. One essen- 
tially counts additional events with a specific signature, here two oppositely charged, almost 
monoenergetic muons, over a smooth background. The cross section for slepton pair pro- 



duction rises as a oc /3 3 , where (3 = y 1 — 4m~ fl /s is the velocity related to the /i# mass. 

The excitation curve as a function of the cms energy, including effects due to QED initial 
state radiation and beamstrahlung, is shown in figure 0. The sensitivity to the width T/x R 
is most pronounced close to the kinematic production limit and diminishes with increasing 
energy. A larger width 'softens' the rise of the cross section with energy. Fits to various mass 
and/or width hypotheses are performed by simulating measurements with a total integrated 
luminosity of 100 fb -1 distributed over 10 equidistant points around s/s = 264 — 274 GeV. 
The data may be collected within a few months of Tesla operation. 
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Figure 2: Observable cross section near threshold of the reaction e R e\ — > Hr^r including 
QED radiation and Beamstrahlung. Curves assume a mass m^ fl = 132.0 GeV and width 
T~ m of 310 MeV (full curve), MeV (dashed curve) and 620 MeV (dashed-dotted curve). 
Measurements correspond to C = 10 fb _1 per point. 

Taking the width from the model prediction r^ fl = 310 MeV, a fit to the threshold curve 
gives a statistical accuracy for the smuon mass of Sm^ R = 90 MeV. This error is considerably 
smaller than the expected width. A two-parameter fit yields = 132.002 1° 130 GeV an< ^ 
^p. R = 311 ^225 MeV. However, both parameters are highly correlated with a correlation 
coefficient of 0.95. Finally, if one may fix the /Ir mass from another measurement, the width 
can be determined to ST^ = ±190 MeV. It should be noted that the scan procedure and 
choice of energy points is by no means optimized. Possibilities to reduce the correlations 
should be studied. 

2.2 Energy spectrum of ^ 

For energies far above threshold, the kinematics of the decay chain of reaction (1) allows to 
identify and to reconstruct the masses of the primary and secondary sparticles. The isotropic 
decays of the scalar muons lead to a flat energy spectrum of the observed final fi in the 
laboratory frame. The endpoints of the energy distribution are related to the masses of the 
jiR and Xi y i a 



In practice the sharp edges of the energy spectrum will be smeared by effects due to detec- 
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tor resolution, selection criteria and in particular initial state radiation and beamstrahlung. 
The results of a simulation at y/s = 320 GeV assuming an integrated luminosity of 160 fb _1 
are shown in figure |3|. One observes a clear signal from p,R pair production above a small 
background of cascade decays X2 ~^ from the reaction e^e\ — > X2XV Contamination 

from chargino or W pair production is completely negligible. 



120 




40 60 80 100 

leptoxi energy Ej [GeV] 



8000 r 
7000 r 
6000 : 
5000 '- 
4000 f- 
3000 : 
2000 '- 
1000 - \- 



20 22 24 26 28 30 32 
lepton Energy E (Gev) 



Figure 3: Left: Energy spectrum of di-muon events from the reaction e. R e\ — > [Ir[Ir and 
the background e^e£ — > X1X2 a ^ \A = 320 GeV assuming £ = 160 fb _1 . Right: Lower 
endpoint region of the /x energy spectrum illustrating the effect of various widths T^ R of 
310 MeV (dots), MeV (dashed) and 620 MeV (dotted) using the tenfold luminosity 

A two-parameter fit to the \x energy spectrum yields masses of m^ R = 132.0 ± 0.3 GeV 
and m x o = 71.9 ±0.2 GeV. The statistical accuracy is of the same size as the expected width 
of the scalar muon. Choosing a different width in the simulation modifies essentially 
the n energy spectrum at the low endpoint and has little impact at higher energies. This is 
illustrated in figure right part, which compares the lower part of the spectrum with the 
predictions of width zero and twice the expected value. The sharp rise is getting smeared 
out with increasing width. With the anticipated luminosity of 160 fbT 1 it may be feasible 
to distinguish these cases. 

2.3 Production of other sparticles 

It should be noted that estimates on the sensitivity of width effects in other slepton produc- 
tion channels can be obtained from the above results by scaling the cross section and taking 
the branching ratios into final states into account. Thus one expects e.g. a gain by a factor 
of ~2 for selectron and sneutrino v e pair production. For the higher mass chargino xt 
and neutralinos ^3,4 mass resolutions of ~ 0.25 — 0.50 GeV may be obtained from threshold 
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scans [f|], where the cross sections rise as a oc (5. The corresponding widths are expected to 
be ~ 2 — 5 GeV (two-body decays in a gauge boson and gaugino) and have certainly to be 
considered. 

3 Conclusions 

The high luminosity of Tesla allows to study the production and decays of the accessible 
Susy particle spectrum. Polarization of both e~ and e + beams is very important to optimize 
the signal and suppress backgrounds. A simulation of slepton production e + e~ — > fi^^Jj 
shows that for precision mass measurements with an accuracy of (9(100 MeV) the widths 
of the primary particles have to be taken into account. Finally, it is worth noting that 
the anticipated mass resolutions from threshold scans or lepton energy spectra can only be 
obtained if beamstrahlung effects are well under control. 
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Radiative Effects on Squark Pair Production at e + e Colliders 

M. DREES, O.J. P. EBOLI, R.M. GODBOLE and S. KRAML 

Abstract 

We study the impact of various radiative effects on the kinematic reconstruction of 
squark pair events at future e + e~ colliders. We focus on the simplest case where both 
squarks decay directly into a stable neutralino, but include both photon radiation off 
the initial state and gluon radiation during the production and/or decay of the squarks. 
These effects change the shapes of the distributions used to determine the squark mass; 
they can therefore also introduce additional systematic uncertainties. 

Among the advantages of e + e~ colliders over hadron colliders are the fairly well denned 
center-of-mass energy in hard (annihilation) events, and the comparative cleanliness of 
the environment. Taken together, these properties make the kinematic reconstruction of 
e + e~ events much easier than that of comparable events at hadron colliders. In particular, 
experiments at e + e~ colliders should be able to measure the masses of new particles (with 
unsuppressed electroweak couplings) with an error of 1% or less Q, either through threshold 
scans, or through fitting kinematic distributions of events well above threshold. The latter 
method is more versatile, since the same data set can well contain several kinds of "new 
physics" events, which can be separated from each other (and from backgrounds) using 
kinematical cuts. Moreover, kinematical reconstruction allows to determine not only the 
masses of the new particles produced in the primary interaction, but also those of their 
decay products. 

However, in order to correctly interpret the information contained in various distributions 
one needs an accurate model of the final state. In particular, if one wants to achieve errors 
of 1% or less, various radiative effects have to be taken into account. These are of special 
importance if the new particles and/or their decay products have strong interactions, since 
a significant fraction of all signal events will then contain hard gluons. Here we study 
the impact of these radiative effects on the measurement of squark masses at e + e~ colliders. 
Squarks are the currently most plausible new particles that have both strong and electroweak 
interactions. 

In the last few years considerable progress has been made [f2| in the accurate calculation of 
total cross sections for squark pair production and of squark branching ratios. In particular, 
one-loop corrections to these quantities from both ordinary QCD and SUSY QCD, as well 
as from Yukawa interactions, are now known. These have been used || to estimate the 
error with which the squark mass can be extracted from a measurement of the total squark 
pair production cross section times branching ratio into a given final state. However, in 
such a "dynamical" determination of the squark mass one has to assume values for all the 
other input parameters that affect the cross section times branching ratio. These include the 
gluino mass and, for third generation squarks, also the masses and mixing angles of squarks, 
charginos, neutralinos, and Higgs bosons 0. 
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In contrast, relatively little attention has been paid to the kinematical determination of 
squark masses at e + e~ colliders. The pioneering work by Feng and Finnell |4[] investigates 
the usefulness of various kinematical distributions, and concludes that experiments at a 
500 GeV collider should be able to determine the mass of 200 GeV squarks with an error 
of < 0.5% using just 20 ftr 1 of data, if all squarks decay directly into a massless quark 
and an invisible (stable or long-lived) neutralino Xv However, their study did not include 
any radiative effects. Here we update their analysis by including initial state radiation of 
photons, as well as the emission of hard gluons during the production (e + e~ — > qq*g) and/or 
decay (q — > qxlg) of the squarks. We also take larger squark masses (rriq ~ 300 GeV) to 
account for recent experimental bounds from the Tevatron , and a correspondingly higher 
center-of-mass energy of 800 GeV. 

For reasons of space we only briefly summarize the ingredients of our analysis; details 
will be given elsewhere. We treat the emission of photons off the initial state in the structure 
function formalism [[J. The differential cross section is then given by 

da = J dx 1 dx 2 f e \ e {x 1 , Vs)f e \ e (x 2 , y/s)da(s = XiX 2 s), (1) 



where a is the cross section in the absence of ISR, and 

f e \ e (x, y/s) = (3 



(2) 



with (3 = (log ^2 — 1 j , is the leading-log resummed effective e ± distribution function. 

Note that we do not include beamstrahlung (which is expected to further smear out the 
peak in the e + e~ luminosity at s — s), since it depends on details of accelerator design. 
Moreover, we conservatively assume that all ISR photons escape detection, even though 
eq.(|2|) is strictly valid only if there are no experimental constraints on the phase space of the 
emitted photons. However, we will see that the main effect of ISR is an overall reduction 
of the cross section by ~ 15%; kinematical distributions are little affected even if all ISR 
photons are invisible. 

We treat the emission of gluons during qq* production as described in 0. In particular, 
we introduce a minimal gluon energy E gyTa \ n to regularize IR divergences. The final results 
will not depend on the value of this regulator after contribution from qq* and qq*g events 
have been added, if -E^min is sufficiently smalLf] In the numerical results presented below we 
take E gMn = 1 GeV. For = 800 GeV and m~ q = 300 GeV this implies that only 18% of 
all squark pairs are produced together with a "hard" gluon. 

The squared matrix element for q — > qx\g can be found in ||. We again have to include 
virtual QCD corrections, in this case to q — > qxi decays, to cancel IR singularities. In this 
case we regularize these singularities by introducing a finite gluon mass m g , which we also 
set to 1 GeV in our numerical examples. For our choice rriq = 300 GeV, = 50 GeV 
this means that nearly 90% of all squark decays produce a "hard" gluon. Altogether more 

*Of course, virtual QCD corrections to qq* production have to be added for this cancellation to work. 
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than 95% of all e + e~ — > qq* — > qqXiXi^ events therefore contain at least one "hard" gluon, 
as defined through our two IR regulators. Note that, unlike for e + e~ — > qq*, virtual QCD 
corrections to q — > qxi decays introduce UV divergences. These cancel only after including 
full SUSY-QCD corrections |§, which depend on the mass of the gluino. We take rrig = 450 
GeV; the shapes of the kinematical distributions we study here are almost independent of 
this choice. 

Due to the emission of hard gluons we have up to 5 visible partons in the final state[]. In 
the numerical results presented below we group these into exactly 2 jets (except for events 
with < 2 partons in the acceptance region defined by |cos(6 l )| < 0.9, which we discard), 
using the kx clustering ("Durham") algorithm ||. We have checked that one obtains very 
similar results if one merges partons with a fixed y cnt parameter, rather than a fixed number 
of final state jets, and only uses the two hardest jets for the kinematical analysis. In order to 
simulate experimental resolutions, we smear the energies (but not directions) of all partons 
before jet merging, with a Gaussian error given by 8(E) = 0.3/vS © 0.01 (E in GeV). 
Finally, we apply a set of cuts in order to suppress backgrounds. We require that the energy 
of each jet exceeds 15 GeV, that the missing exceeds 56 GeV (taken from [ffl after scaling 
up from -y/i = 500 GeV to 800 GeV), and that the acoplanarity angle between the two jets 
exceed 30°. 

In Figs. la,b we show the resulting distribution in the jet energy (la) and in the variable 
^g.min (lb) defined below. In the absence of cuts, radiative effects and energy smearing, 
the jet energy distribution should be constant (flat) between two kinematical endpoints, 
and zero elsewhere. The cuts distort this simple shape, producing a small peak near the 
lower edge of the spectrum; this is shown by the solid curve in Fig. la. The histograms 
show "experimental" distributions, based on 3,500 events before cuts; for a single q = 2/3 
SU(2)— singlet (ur) squark, this would require an integrated luminosity of about 205 fb"" 1 
(including ISR and QCD corrections, and assuming all squarks decay directly into qxT)] the 
total cross section is 17.0 (11.7) fb before (after) cuts.| The dotted histogram still does 
not include any radiative effects, but includes a finite energy resolution as described above; 
clearly this has little effect on the shape of the spectrum. As mentioned earlier, including 
ISR (dashed histogram) reduces the cross section, but again does not distort the shape of 
the spectrum very much; in particular, there is only little "leakage" beyond the nominal 
endpoints. In contrast, including gluon emission (solid histogram) does change the shape of 
the distribution. In particular, there are now quite a few entries below the lower nominal 
endpoint. In most of these events one of the quarks falls out of the acceptance region, so 
that one of the jets is entirely made up of gluons. QCD corrections also increase the total 

^We allow gluon emission during the production and both decays simultaneously, i.e. we include events 
with 2 or 3 gluons. These contributions are formally of NNLO. However, since production and decays are 
independent processes, up to terms 0(Tq/y/s), other (as yet unknown) NNLO contributions cannot cancel 
these known contributions. 

■'■The total cross section for ul pair production for the same mass is 26.8 fb. However, in most models 
one expects SU(2) doublet squarks to predominantly decay into charginos and heavier neutralinos, rather 
than directly into the LSP. 
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Figure 1: The jet energy (left) and 

^g,min (right) distributions after cuts. The solid curves show 
the "ideal" distributions in the absence of energy smearing and radiative effects. In the dotted, 
dashed and solid histograms, resolution smearing, ISR, and gluon radiation are "switched on" 
successively. See the text for further details. 

The variable min plotted in Fig. lb has been introduced in Q . It is the minimal possible 
squark mass, if one assumes a fixed value of (no ISR), and if the final state consists of 
exactly two massless quarks and two neutralinos with equal (and known) mass.[] None of 
these assumptions hold in our case. In fact, in some cases the reconstruction described in j| 
is impossible, i.e. some trigonometrical function acquires a value exceeding unity; we discard 
these events. We nevertheless find that radiative effects introduce only a modest amount of 
"leakage" beyond the nominal endpoint of the distribution, which is at rriq. Note, however, 
that these effects significantly broaden the distribution, i.e. many events have migrated to 
lower values of m^ mill (compare the solid and dotted histograms). Since events in the peak 
contribute most to the determination of rriq [[| , this broadening is expected to increase the 
error on rriq. 

Fig. 2 shows that a fit of the m,- im in distribution nevertheless yields a smaller statistical 
error for rriq than a fit of the Ej ct distribution does. In this figure we compare a mock data 
set for rriq = 300 GeV, based on an integrated luminosity of just 50 fb" 1 (852 events before 
cuts), with "template" distributions, which have been computed for 13 different values of 
rriq, leaving all other input parameters unchanged, and applying the cuts described above. 

§Note that QCD corrections markedly increase the acceptance of the cuts. It is therefore also important 
to include their effect on the kinematics when trying to extract rriq from measurements of the qq* production 
cross section. 

^The neutralino mass is expected to be known from analyses of data at lower energy, e.g. from chargino 
pair events. 
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This gives 13 different values of x 2 '-, the open (filled) squares have been computed from 
the E ]et (mg )m i n ) distributions. A parabola is then fitted to x 2 { m q)- The minimum of this 
parabola gives the "measured" value m^o of rriq, while the (la) error is computed from 
X 2 ( m q,o =t 8 m q) = Xmin + 1- Note that the \ 2 values are only based on the shapes of the 
distributions, i.e. the "data" have been normalized to the template before computing \ 2 
for a given assumed value of rriq. Fig. 2 therefore shows the results of purely kinematical 
determinations of rriq. 




Figure 2: The x 2 values computed by comparing a mock data set with 13 different template 
distributions. The "data" are for rriq = 300 GeV, while the templates are for values of rriq between 
288 and 312 GeV. The open (full) boxes are for the -Ej e t (w^mm) distributions, and the curves are 
parabolic fits to the corresponding functions x 2 { m q)- These parabolas determine the "measured" 
rriq, including its error. 



We see that even after radiative effects have been included, a fit of the m,- m i n distribution 
determines the squark mass with a statistical error of well under 1% already with 50 fb~ x 
of data, which corresponds to only about 1 month of running time for the planned TESLA 



collider |]10fl . However, radiative effects also introduce possible new systematic errors, e.g. 
due to the choice of scale in a s . (In our calculation we used for the corrections to qq* 
production, and rriq in q decays.) Moreover, we haven't included any hadronization effects 
yet. Note that the lighter scalar top eigenstate t\, which is likely to be the lightest, and 
hence most easily accessible squark, might well hadronise before it decays [|7|. In addition, 
the massless quarks and gluons in the final state will hadronise into jets with finite masses. 
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Finally, the error on the assumed LSP mass also propagates into the error on rriq H. We 
plan to investigate these issues in a future publication. 
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Spin correlations and phases for SUSY particle searches at e + e 

colliders 

N. GHODBANE 

1 Introduction 

Since one expects high luminosities for the next generation of linear colliders (e.g. ~ 500/6 -1 
for the TESLA project), one can use beam polarization to reduce the standard model back- 
grounds and use the polarization dependence of the cross sections to study specific SUSY 
parameters. Moreover, as it has been stressed by several authors [l|], spin correlations play 
a major role in the kinematic distributions of final particles. For this reason, we upgraded 
the SUSYGEN Monte Carlo generator^ [f| . 

Here, we firstly study the spin correlations effects in the gaugino and the stau searches. 
Then, in a second part, we show how the newly reconsidered CP violating phases arising 
from MSSM can affect the chargino searches. 

2 Beam polarization and spin correlations 

As it has been stressed by several papers |l| the study of the angular distribution of e 
produced in e + e _ — > xlXi ~^ XiXi e+e ~ will give valuable information concerning the neu- 
tralino nature and then enable MSSM parameter extraction. Here we consider the two 
scenarios (A: M 2 = 78 GeV/c 2 , fi = -250 GeV/c 2 , tan/3 = 2) and (B: M 2 = 210 GeV/c 2 , 
fj, = -60 GeV/c 2 , tan/3 = 2). In model A, xl is Bino like (m^o = 42.5 GeV/c 2 ,x? = 
+0.985 + 0.171U 3 - 0.09^° + 0.04#2°) whereas xl is Wino like (m^o = 91.9 GeV /c 2 , xl = 
+0.145 - 0.95IU 3 - 0.28F? + 0.05F 2 °). In model B, the two lightest neutralinos are Hig- 
gsino like (m^o = 55.1 GeV/c 2 ,x? = +0.165 - 0.09IU 3 - 0.78/7? - 0.60# 2 °) and (m^o = 
88.9 GeV /c 2 , xl = +0.21B - 0.24IU 3 + 0.625° - 0.71iJ 2 ). For each one of these two mod- 
els, we considered the common sfermions mass at the GUT scale, mo, being equal to 80 
GeV/c 2 and 200 GeV/c 2 . 

The right side of figure [I] illustrates the effect of spin-correlations in angular distributions 
of e~, decay product of xl- One can notice that the angular distribution of the final leptons 
depends strongly on the sfermion mass parameter m . For small m of inclusion of spin 
correlations gives an effect ~ 20% in the angular distributions. For higher selectron masses, 
the total cross section is smaller, but the spin correlation effects appear to be more important 
~ 30%. If the two neutralinos are Higgsino like, the effects are negligible (~ 0.7%) (see figure 
0). Moreover these angular distributions depend strongly on the center of mass energy. 

Recent studies have shown that r polarization effects yield valuable information for 
the MSSM parameters e.g for tan (3, the nature of the LSP and the mixing angle Of. 

*A description of the Monte Carlo generator can be found in the Higgs working group report. 
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Figure 1: The figure on the left side shows the evolution of the cross section associated to the 
production of X1X2 f° r different beam polarizations. The figure on the right side shows the 
da jd cos 9 distribution for the final e~ decay product of the X2 with (solid line) and without 
spin correlations (dashed line) for two assumptions on the mo parameter (scenario A). 
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Figure 2: The figure on the left side shows the da/d cos 9 distribution for the final e~ 
decay product of the X2 with (solid line) and without spin correlations (dashed line) for two 
assumptions on the mo parameter assuming the neutralino being bino like (scenario A). The 
Higgsino like case (scenario B and figure on the right) shows no strong dependence on the 
spin correlations. 



Figure |3| shows the momentum distribution of pions produced by e + e~ — > f^fi — > 
T+ Xi T Xi n+ ^TXi n ~^TXv The distributions have been plotted assuming two scenarios 
for the stau mixing angle 0^ = and ir, and Xi nature (Bino, Higgsino). One can see that 
the final particle distributions will give access to the tau polarization V T , to tan/5, and to 
Xi nature and through them, to the MSSM parameters. 
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Figure 3: Momentum distributions of n, decay product of r produced in the process e + e — » 
f^f^ — > T + x'j ) T _ Xi- The distributions have been plotted for two hypotheses concerning the 



stau mixing angle Of (tl and tr) and for two hypotheses concerning the neutralino Xi nature. 



3 Phases in supersymmetry searches 

In the MSSM, there are new potential sources of CP non conservation Complex CP 
violating phases can arise from several parameters present in the MSSM Lagrangian: the 
Higgs mixing mass parameter /i, the gauginos masses Mj, the trilinear couplings A^. Exper- 
imental constraints on these CP violating phases come from the electric dipole moment of 
the electron and the neutron ||. 

Figure ^ (left side) shows the chargino pair production cross section variation in terms 
of the phase associated to the \x parameter, for several values of the sneutrino mass 
and for a value of tan/5 = 1.5. One sees that there is a local minimum between the two 
extreme values ($ M = 0, it), which are tested at LEP searches (fi > 0, fi < 0), but also that it 
is not so deep as to raise doubts on the exhaustiveness of "phaseless" searches. Further the 
electric dipole moment experimental upper limit {E e * v < 4.3 • 10 _27 e ■ cm) will constrain 
these phases (figure on the right side) and rule out many scenarios for which the smallest 
cross section for chargino pair production is obtained for a M parameter different from 
and 7r. 

At the linear collider, one cannot neglect this strong dependence of the cross section on 
phases. 
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The three-leptons signature from resonant sneutrino production 

at the LHC 

G. MOREAU, E. PEREZ and G. POLESELLO 

Abstract 

The resonant production of sneutrinos at the LHC via the R-parity violating cou- 
plings \' i j k LiQjD'j. is studied through its three-leptons signature. A detailed particle 
level study of signal and background is performed using a fast simulation of the ATLAS 
detector. Through the full reconstruction of the cascade decay, a model- independent 
and precise measurement of the masses of the involved sparticles can be performed. 
Besides, this signature can be detected in a large part of the SUSY parameter space 
and for wide ranges of values of several A' ijfc coupling constants. 

1 Introduction 

In extensions of the Minimal Supersymmetric Standard Model (MSSM) where the so-called 
R-parity symmetry is violated, the superpotential contains some additional trilinear cou- 
plings which offer the opportunity to singly produce supersymmetric (SUSY) particles as 
resonances. The analysis of resonant SUSY particle production allows an easier determina- 
tion of the these R-parity violating (R p ) couplings than the displaced vertex analysis for the 
Lightest Supersymmetric Particle (LSP) decay, which is difficult experimentally especially 
at hadronic colliders. 

In this paper, we study the sensitivity provided by the ATLAS detector at the LHC 
on singly produced charginos via the A 211 coupling, the main contribution coming from the 
resonant process pp — > — > xfl^ T - At hadron colliders, due to the continuous energy distri- 
bution of the colliding partons, the resonance can be probed over a wide mass range. We have 
chosen to concentrate on X'^LiQjD^ interactions since Xl^UfD^DI couplings lead to mul- 
tijet final states with large QCD background. Besides, we focus on A' 211 since it corresponds 
to first generation quarks for the colliding partons and it is not severely constrained by low 
energy experiments: A 211 < 0.09 (for m = 100 GeV) [JJ. We consider the cascade decay lead- 
ing to the three-leptons signature, namely Xi ~^ XilpVp (with l p = e, fi), Xi ~^ l^ud, Jiud. 
The main motivation lies in the low Standard Model background for this three-leptons final 
state. The considered branching ratios are typically of order B(xt Xi^p^p) ~ 22% (for 
mi,mg,m^o > rn~±) and B(xi ~^ fiud) ~ 40%— ~ 70%. 

2 Mass reconstruction 

The clean final state, with only two hadronic jets, three leptons and a neutrino allows the 
reconstruction of the v decay chain and the measurement of the x?> Xi an d ^ masses. We 
perform full analysis for the following point of the MSSM: Mi = 75 GeV, M 2 = 150 GeV, 



38 



fi = -200 GeV, tan/? = 1.5, A t = A b = A T = 0, m~ f = 300 GeV and for A' 211 =0.09. For this 
set of MSSM parameters, the mass spectrum is: m^o = 79.9 GeV, m-± = 162.3 GeV and 
the total cross-section for the three-leptons production is 3.1 pb, corresponding to ~ 100000 
events for the standard integrated luminosity of 30 fb _1 expected within the first three years 
of LHC data taking. 

The single chargino production has been calculated analytically and implemented in a 
version of the SUSYGEN MonteCarlo || modified to include the generation of pp processes. 
The generated signal events were processed through the program ATLFAST [Q, a parame- 
terized simulation of the ATLAS detector response. 

First, we impose the following loose selection cuts in order to select the considered final 
state and to reduce the Standard Model (SM) background (see Section p.l|) : (a) Exactly 
three isolated leptons with p\ > 20 GeV, p^ 3 > 10 GeV and |^| < 2.5, (b) At least two of 
the three leptons must be muons, (c) Exactly two jets with p? > 15 GeV, (d) The invariant 
mass of any pair must lie outside ±6.5 GeV of the Z mass. 

The three leptons come in the following flavor-sign configurations (+ charge conjugates): 
(1) fi~e + fi + (2) /i~e + /i~ (3) fi~fi + fi + (4) where the first lepton comes from the 

z/^, the second one from the W and the third one from the Xi decay. As a starting point 
for the analysis, we focus on configuration (1) where the muon produced in the x? decay is 
unambiguously identified as the one with the same sign as the electron. The distribution of 
the //-jet-jet invariant mass exhibits a clear peak over a combinatorial background, shown 
on the left side of Figure [TJ. After combinatorial background subtraction (right of Figure [I]) 
an approximately Gaussian peak is left, from which the x? mass can be measured with a 
statistical error of ~ 100 MeV. The combinatorial background is due to events where one 
jet from \\ decay is lost and a jet from initial state radiation is used in the combination, 
and its importance is reduced for heavier sneutrinos or neutralinos. Once the position of the 
Xi mass peak is known, the reconstructed x? statistics can be increased by also considering 
signatures (2), (3) and (4), and by choosing as the x? candidate the muon-jet-jet combination 
which gives invariant mass nearest to the peak measured previously using events sample (1). 
For further reconstruction, we define as x? candidates the /^-jet-jet combinations with an 
invariant mass within 12 GeV of the measured x? peak, yielding a total statistics of 6750 
events for signatures (1) to (4) for an integrated luminosity of 30 fb _1 . 

For Xi reconstruction we consider only configurations (1) and (2), for which the charged 
lepton from W decay is unambiguously identified as the electron. The longitudinal momen- 
tum of the neutrino from the W decay is calculated from the missing transverse momentum 
of the event (p^) and by constraining the electron-neutrino invariant mass to the W mass. 
The resulting neutrino longitudinal momentum has a twofold ambiguity. We therefore build 
the invariant W — x? mass candidate using both solutions for the W boson momentum. 
The observed peak, represented on the left side of Figure [| can be fitted with a Gaussian 
shape with a width of ~ 6 GeV. Only the solution yielding the xf mass nearer to the mea- 
sured mass peak is retained, and the x^ candidates are defined as the combinations with an 
invariant mass within 15 GeV of the peak, corresponding to a statistics of 2700 events. 

Finally, the sneutrino mass is reconstructed by taking the invariant mass of the Xi 
candidate and the leftover muon (Figure ^, right). The v mass peak has a width of ~ 10 GeV 
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and 2550 events are counted within 25 GeV of the measured peak. 




Figure 1: fi-jet-jet invariant mass for events in configuration (1) (see text) before (left) and 
after (right) background subtraction. 



3 Analysis reach 

3.1 Standard Model background 

We consider the following SM processes for the evaluation of the background to the three- 
leptons signature: (1) it production, followed by t — > Wb, where the two W and one of the 
b quarks decay leptonically, (2) WZ production, where both bosons decay leptonically, (3) 
Wt production, (4) Wbb production, (5) Zb production. These backgrounds were generated 
with the PYTHIA Monte Carlo |4[] , and the ONETOP parton level generator || , and passed 
through the ATLFAST package §. 

We apply to the background events the loose selection cuts described in Section 0, and 
in addition we reject the three same-sign muons configurations which are never generated 
by our signal. The background to the sneutrino decay signal is calculated by considering 
the events with a /x-jet-jet invariant mass in an interval of ±15 GeV around the x? peak 
measured for the signal. In order to optimize the signal to background ratio only events 
containing three muons (configurations (3) and (4)), which are less likely in the Standard 
Model, are considered. In each event two combinations, corresponding to the two same-sign 
muons, can be used for the Xi reconstruction. Both configurations are used when counting 
the number of events in the peak. In most cases, however, the difference in mass between 
the two combinations is such that they do not appear in the same peak region. 
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Figure 2: Left: Invariant mass of the W with the x\ candidate. Right: Invariant mass of 
the third lepton in the event with the xf candidate. 



3.2 Supersymmetric background 

The pair production of SUSY particles through standard i? p -conserving processes represents 
another source of background. A study based on the HERWIG 6.0 MonteCarlo H has 



shown that all the SUSY events surviving the cuts described in Section |3.1| are mainly 
from pp — > x + X reactions (x being either a chargino or a neutralino and X any other 
SUSY particle), and that the SUSY background decreases as the x^ an d X° masses increase. 
This behavior is due to the combination of two effects: the x + X production cross-section 
decreases with increasing x mass, and the probability of losing two of the four jets from the 
decays of the two Xi m the event becomes smaller as the x and x? masses increase. The 
SUSY background is only significant for xf masses lower than 200 GeV. 

Besides, it can be assumed that the Xi mass will be derived from inclusive Xi reconstruc- 
tion in SUSY pair production as shown in |7|] and ||. Hence, even in the cases where a 
significant Xi peak can not be observed above the SUSY background, we can proceed to the 
further steps in the kinematic reconstruction. The strong kinematic constraint obtained by 
requiring both the correct x\ mass and a peak structure in the x? — W invariant mass will 
then allow to separate the single sneutrino production from other SUSY processes. 

Therefore, only the Standard Model background is considered in the evaluation of the 
analysis reach presented below. 
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3.3 Reach in the mSUGRA parameter space 



In Figure |3|, we show the regions of the mo-mi/2 plane where the signal significance exceeds 
5 a (-7= > 5 with S = Signal and B = SM Background) after the set of cuts described 
in Section |3.1| has been applied, within the mSUGRA model. The full mass reconstruction 
analysis of Section ||] is possible only above the dashed line parallel to the m axis. Below 
this line the decay Xi ~ * X^W ± * s kinematically closed, and the W mass constraint can not 
be applied to reconstruct the neutrino longitudinal momentum. 




Figure 3: 5a reach in the mo — mi/2 plane for Aq = 0, tan/3 = 2, \i < (left) and fi > 
(right) and three different choices of the A' 211 coupling, with an integrated luminosity of 30 fb^ 1 
at the LHC. A signal of at least ten events is required. The hatched region at the upper left 
corresponds to m„ < rn^t. The cross-hatched region for low m 1( / 2 gives the kinematical limit 

for the discovery of xf or I by LEP running at y/s = 196 GeV j^J. The dotted line shows 
the region below which the xf decays to a virtual W. 



The basic feature in Figure |3] is a decrease of the sensitivity on A' 211 as m increases. This 
is due to a decrease of the partonic luminosity as increases. The sensitivity on A' 211 is 
also observed to decrease as m~± approaches m^. There are two reasons. First, in this region 
the phase space factor of the decay v — > xt nT following the resonant sneutrino production 
is suppressed, thus reducing the branching fraction. Secondly, as the and the xt become 
nearly degenerate the muon from the decay becomes on average softer, and its pr can fall 
below the analysis requirements. In the region m-± > mp, shown as a hatched region in the 
upper left of the plots, the resonant sneutrino production contribution vanishes and there is 
essentially no sensitivity to A' 211 . Finally, the the sensitivity vanishes for low values of mi/2- 
This region, below the LEP 200 kinematic limit for Xi detection, corresponds to low values 
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of the Xi mass. In this situation the two jets from the Xi decay are soft, and one of them 
is often below the transverse momentum requirement, or they are reconstructed as a single 
jet. 

For high tan /3, the three-lepton signature is still present, but it may be produced through 
the decay chain xf ~* T\v T , followed by f\ — ► r%5- The full kinematic reconstruction becomes 
very difficult, but the signal efficiency is essentially unaffected, as long as the mass difference 
between the lightest f and the Xi is larger than ~ 50 GeV. For a smaller mass difference the 
charged lepton coming from the r decay is often rejected by the analysis cuts. 

4 Conclusion 

In conclusion we have shown that if minimal supersymmetry with R-parity violation is 
realized in Nature, the three-leptons signature from resonant sneutrino production will be a 
privileged channel for precisely measuring sparticle masses in a model-independent way as 
well as for testing a broad region of the mSUGRA parameter space. 

This signature can lead to a high sensitivity to the A' 211 coupling and should also allow to 
probe an unexplored range of values for many other R p couplings of the type A' ljfe and X' 2 j k . 
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Resonant slepton production at the LHC 

H. DREINERQ, P. RICHARDSONQ, and M. H. SEYMOURf| 



Abstract 

We consider the resonant production of charged sleptons at the LHC via R-parity 
violation (R p ) followed by gauge decays to a charged lepton and a neutralino which 
then decays via R p . This gives a signature of two like-sign charged leptons. In the 
simulation we include the full hadronization via Monte Carlo programs. We find a 
background, after cuts, of 5.1 ± 2.5 events for an integrated luminosity of 10/6 -1 . A 
preliminary study of the signal suggests that couplings of 2 x 10 -3 for a smuon mass 
of 223GeV and smuon masses of up to 540GeV for couplings of 10 -2 can be probed. 

1 Introduction 

In R-parity violating (R p ) models the single resonant production of charged sleptons in 
hadron-hadron collisions is possible. The most promising channels for the discovery of these 
processes, at least with small R p couplings, involve the gauge decays of these resonant 
sleptons. In particular if we consider the production of a charged slepton, this can then 
decay to give a neutralino and a charged lepton, i.e. the process 

u + d — >£+ — >£ + + x°- (1) 

In addition to this s-channel process there are t-channel processes involving squark exchange. 
The neutralino decays via the crossed process to give a charged lepton, which due to the 
Majorana nature of the neutralino can have the same charge as the lepton from the slepton 
decay. We therefore have a like-sign dilepton signature which we expect to have a low 
Standard Model background. 

2 Backgrounds 

The dominant Standard Model backgrounds to this process come from 

• Gauge boson pair production, i.e. production of ZZ or WZ followed by leptonic decays 
of the gauge bosons with some of the leptons not being detected. 

• tt production. Either the t or t decays semi-leptonically, giving one lepton. The second 
top decays hadronically. A second lepton with the same charge can be produced in a 

*E-mail address: drciner@v2.rl.ac.uk 

^E-mail address: p.richardsonl@physics.ox.ac.uk 

•^E-mail address: M.Seymour@rl.ac.uk 
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semi-leptonic decay of the bottom hadron formed in the hadronic decay of the second 
top, i.e. 



t -> W + b^e+z7 e b, 

t -> W"b -> qqb, b -> e + z7 e c. (2) 

• bb production. If either of these quarks hadronizes to form a s meson this can mix 
to give aBjj s . This means that if both the bottom hadrons decay semi-leptonically the 
leptons will have the same charge as they are both coming from either b or b decays. 

• Single top production. A single top quark can be produced together with a b quark 
by either an s- or t-channel W exchange. This can then give one charged lepton from 
the top decay, and a second lepton with the same charge from the decay of the meson 
formed after the b quark hadronizes. 

• Non-physics backgrounds. There are two major sources: (i) from misidentifying the 
charge of a lepton, e.g. in Drell-Yan production, and (ii) from incorrectly identifying 
an isolated hadron as a lepton. This means that there is a major source of background 
from W production with an additional jet faking a lepton. 

Early studies of like-sign dileptons at the LHC [[IJ] only studied the backgrounds from 
heavy quark production. It was found that by imposing cuts on the transverse momentum 
and isolation of the leptons the heavy quark backgrounds could be significantly reduced. 
However more recent studies of the like-sign dilepton production at the LHC and the 
Tevatron ||,[|| suggest that a major source of background to like-sign dilepton production 
is from gauge boson pair production and from fake leptons. Here we will consider the 
backgrounds from gauge boson pair production as well as heavy quark production. The 
study of the non-physics backgrounds (e.g. fake leptons) requires a full simulation of the 
detector and it is therefore beyond the scope of our study. In particular the background 
from fake leptons cannot be reliably calculated from Monte Carlo simulations and must 
be extracted from data ||. We can use the differences between the Jl p signature we are 
considering and the MSSM signatures considered in [^J to reduced the background from 
gauge boson pair production. 

We impose the following cuts 

• A cut on the transverse momentum of the like-sign leptons px > 40 GeV. 

• An isolation cut on the like-sign leptons so that the transverse energy in a cone of 
radius R = \J A<p 2 + Arf = 0.4 about the direction of each lepton is less than 5 GeV. 

• A cut on the transverse mass, M% = 2|p T J |p T J(l— cos where px e is the transverse 
momentum of the charged lepton, p? v is the transverse momentum of the neutrino, 
assumed to be all the missing transverse momentum in the event, and A<p£ U is the 
azimuthal angle between the lepton and the neutrino, i.e. the missing momentum in 
the event. We cut out the region where 60GeV < Mt < 85GeV. 
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Figure 1: Effect of the isolation cuts on the tt and bb backgrounds. The dashed line gives 
the background before any cuts, the solid line shows the effect of the isolation cut described 
in the text. The dot-dash line gives the effect of all the cuts. 



• A veto on the presence of a lepton in the event with the same flavor but opposite charge 
(OSSF) as either of the leptons in the like-sign pair if the lepton has pr > 10 GeV and 
which passes the same isolation cut as the like-sign leptons. 

• A cut on the missing transverse energy, E^ iss < 20 GeV . 

While these cuts were chosen to reduce the background we have not attempted to op- 
timize them. The first two cuts are designed to reduce the background from heavy quark 
production. As can be seen in Fig. [I], these cuts reduce this background by several orders of 
magnitude. The remaining cuts are designed to reduce the background from gauge boson 
pair, in particular WZ, production which is the major source of background after the imposi- 
tion of the isolation and px cuts. The transverse mass cut is designed to remove events with 
leptonic W decays as can be seen in Fig. 0a. The veto on the presence of OSSF leptons is 
designed to remove events where one lepton from the dilepton pair comes from the leptonic 
decay of a Z boson. The missing transverse energy cut again removes events with leptonic 
W decays, this is mainly to reduce the background from WZ production, as seen in Fig. 0b. 
The effect of these cuts on the heavy quark and gauge boson pair backgrounds are shown in 
Figs. |] and 0, respectively. 

The backgrounds from the various processes are summarized in Table [|. The simula- 
tions of the bb, tt and single top production were performed using HERWIG6.1 [[J. The 
simulations of gauge boson pair production used PYTHIA6.1 ||. The major contribution to 
the background comes from WZ production; the major contribution to the error comes from 
bb. For the bb simulation we have required a parton-level cut of 40 GeV on the transverse 
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Figure 2: Transverse mass and missing transverse energy in WZ events 
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Figure 3: Effect of the isolation cuts on the WZ and ZZ backgrounds. The dashed line gives 
the background before any cuts, the solid line shows the effect of the isolation cut described 
in the text. The dot-dash line gives the effect of all the cuts. 
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Background Process 


Number of Events 


After pt cut 


After isolation and pr cuts 


After all cuts 


WW 


2.8 ±0.6 


0.0 ±0.1 


0.0 ±0.1 


wz 


226 ±3 


189 ±3 


4.1 ±0.5 


zz 


50.4 ±0.9 


40.6 ±0.8 


0.9 ±0.1 


tt 


(4.8 ±0.3) x 10 3 


0.34 ±0.14 


0.06 ±0.06 


bb 


(5.69 ±0.8) x 10 4 


0.0 ±2.4 


0.0 ±2.4 


Single Top 


11.5 ±0.3 


0.0 ±0.008 


0.0 ±0.008 


Total 


(6.2 ±0.8) x 10 4 


230 ±4 


5.1 ±2.5 



Table 1: Backgrounds to like-sign dilepton production at the LHC. The numbers of events 
are based on an integrated luminosity of 10 fb _1 . We used the cross sections from the Monte 
Carlo simulation for bb and single top production, the next-to-leading order cross section for 
gauge boson pair production from and the next-to-leading order with next-to-leading-log 
resummation cross section from || for tt production. We estimate an error on the cross 
section from the effect of varying the scale between half and twice the hard scale, apart from 
gauge boson pair production where we do not have this information for the next-to-leading 
order cross section. The error on the number of events is then the error in the cross section 
and the statistical error from the simulation added in quadrature. 

momentum of the bottom quarks. This should not affect the results provided we impose a 
cut of at least 40 GeV on the px of the leptons. We also forced the B meson produced to 
decay semi-leptonically. In events where there was one B^ s meson this meson was forced to 
mix, if there was more than one B° d s then one of the mesons was forced to mix and the others 
forced to not mix. Even with these cuts it is impossible to simulate the full luminosity with 
the resources available, due to the large cross section for bb production. This gives the large 
error on the estimate of this background. 

3 Signal 

We used HERWIG6.1 to simulate the signal. This version includes the resonant slepton 
production, including the t-channel diagrams, and the R-parity violating decay of the neu- 
tralino including a matrix element for the decay |J. We will only consider first generation 
quarks as the cross sections for processes with higher generation quarks are suppressed by 
the parton distributions. There are upper bounds on the R p couplings from low energy 
experiments. The bound on A'm from neutrino-less double beta decay [|ITJ,|TTJ is very strict 
so we consider muon production via the coupling A^n, which has a much weaker bound, 

A - < 059 x (issb) • < 3) 

from the ratio R v = T(n — > eu)/T(7i — > fiu) 0,0. 



18 



a. Number of Signal Events b. Efficiency of Cuts 
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Figure 4: Number of signal events passing the cuts and the efficiency for Mi/ 2 = 300GeV, 
A = 300GeV, tan/3 = 2, sgn/i = +, with the Rp coupling A' 2 n = 0.01. The dashed line 
gives the number of events needed for a 5<r discovery. 



We have performed a scan in Mq using HERWIG with the following SUGRA parameters, 
Mi /2 = 300GeV, A = 300GeV, tan (3 = 2, sgn/x = +, and with the R p coupling X' 211 = 0.01. 
The number of events which pass the cuts given in Section ^| are shown in Fig. |]a, while the 
efficiency of the cuts, i.e. the fraction of the signal events which have a like-sign dilepton 
pair passing the cuts, is shown in Fig.|]b. The dip in the efficiency between 140GeV < 
Mo < 180GeV is due to the resonant production of the second lightest neutralino becoming 
accessible. Just above threshold the efficiency for this channel is low due to the low pr of 
the lepton produced in the slepton decay. 

If we conservatively take a background of 7.6 events, i.e. la above the central value of 
our calculation, a 5<r fluctuation of the background would correspond to 20 events, using 
Poisson statistics. This is given as a dashed line in Fig. |]a. As can be seen for a large range 
of values of Mq resonant slepton production can be discovered at the LHC, for A 211 = 0.01. 
The production cross section depends quadratically on the R p Yukawa coupling and hence 
it should be possible to probe much smaller couplings for small values of M . 

As can be seen in Fig. |5], at this SUGRA point the sdown mass varies between 622GeV at 
M = 50GeV and 784GeV at M = 500GeV. The corresponding limit on the coupling A^n 
varies between 0.37 and 0.46. We can probe couplings of A^n = 2 x 10~ 3 for Mo = 50GeV 
which corresponds to a smuon mass of 223GeV, and at couplings of A^n = 10~ 2 we can probe 
values of Mo up to 500GeV, i.e. a smuon mass of 540GeV. This is more than an order of 
magnitude smaller than the current upper bounds on the R p coupling given above for these 
values of M . This is a greater range of couplings and smuon masses than can be probed at 
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Figure 5: Masses of the left smuon, solid line, and the right sdown, dashed line, as a function 
of M for M V 2 = 300GeV, A = 300GeV, tan (3 = 2, sgn/x = +. 



the Tevatron |13|]. The backgrounds are higher at the LHC but this is compensated by the 
higher energy and luminosity leading to significantly more signal events. 



4 Conclusions 

We have considered the backgrounds to like-sign dilepton production at the LHC and find 
a background after cuts of 5.1 ± 2.5 events for an integrated luminosity of 10/6 -1 . This 
means, taking a conservative estimate of the background of 7.6 events, that 20 events would 
correspond to a 5cr discovery. For a full analysis however, non-physics backgrounds must 
also be considered. 

A preliminary study of the signal suggests that an efficiency for detecting the signal in 
excess of 20% can be achieved over a range of points in SUGRA parameter space. At the 
SUGRA point studied this means we can probe R p couplings of 2 x 10~ 3 for a smuon mass 
of 223GeV and up to smuon masses of 540GeV for couplings of 10 -2 , and higher masses for 
larger couplings. 

A more detailed scan of SUGRA parameter space for this signal remains to be performed. 
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Xi reconstruction in mSUGRA models with R— parity breaking 

LQD term 

L. MEGNER and G. POLESELLO 
Abstract 

The reconstruction of the Xi LSP through the decay to lepton-jet-jet is studied 
in the framework of the mSUGRA model. A detailed particle level analysis is per- 
formed on events generated with the HERWIG 6.0 Monte Carlo and passed through 
the ATLAS fast simulation program ATLFAST. The extraction of the Xi mass peak 
and techniques for subtracting the combinatorial background are demonstrated on four 
example points in the mSUGRA parameter space. 



1 Introduction 

The R— parity violating ( R p ) extension of the MSSM contains the following additional terms 
in the superpotential, which are trilinear in the quark and lepton superfields, 

W R . odd = ]T (\KjkLiLjE* + KjkLiQjDt + ^! jk U'Dpt) , (1) 

where i, j, k are flavor indices. The presence of R- violating terms will mainly manifest itself 
in two ways: by production of single sparticles, or by the fact that the LSP is not stable and 
decays to standard model particles. See for a recent review on the subject. 

The reconstruction of the x? from its decay products in SUSY particle pair production 
offers the possibility of detecting R— parity violation for values of the .R-violating couplings 
too small to yield a significant cross-section for single sparticle production. The x? decay can 
be studied for all three coupling types, A, A' and A" || ||. The hypothesis of a dominating 
A' is the most favorable case for the kinematic reconstruction of SUSY cascade decays, 
if the decay of the x? into a charged lepton has a sizeable branching fraction. In this 
case, the charged lepton can be used as the initiator of the mass reconstruction, yielding 
a reduced combinatorial background with respect to three-jet case, and all three particles 
from the x° decay are reconstructed in the detector, allowing to fully reconstruct the x? mass 
peak. In a pioneering work [|J it was explicitly demonstrated that with appropriate selection 
criteria on jets and leptons a clear x\ peak could be observed in the lepton-jet-jet invariant 
mass over an acceptable combinatorial background. That work was based on x? decaying 
with 100% branching fraction to eqq', with the momenta of the decay products distributed 
according to the phase space. The present study is based on a detailed implementation of 
the R— parity violating decays in the HERWIG 6.0 Monte Carlo including matrix elements 
for the Xi decay [7j . We aim to to verify if the x? peak can be observed for a broad range of 
SUSY models, and to develop techniques for evaluating and subtracting the combinatorial 
background, as a starting point for performing precision measurements of the parameters of 
the underlying model. 
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2 The SUSY model 



We work in the framework of the minimal Supergravity inspired model (mSUGRA) as im- 
plemented in the ISASUSY Monte Carlo [||]. The SUSY events were generated with version 

6.0 of the HERWIG Monte Carlo 0, which implements all the R— parity conserving two-to- 
two SUSY processes. This version of HERWIG also includes the simulation of all R— parity 
violating sparticle decays H. The produced events were then passed through ATLFAST, a 
fast simulation of the ATLAS detector || . 

We assumed for this study a single dominant R-violating coupling A / 111 =0.01. The results 
obtained are independent of the precise value of the A' coupling, as long as it induces a prompt 
Xi decay. The results are valid for all the couplings of the type A^ fe with i = 1,2. They 
should be taken with some care if j or k is equal to three, since in this case a b quark is 
among the Xi decay products. 

In models with a single dominant A' term Xi has two decay modes: Xi ~^ I QQ' an d 
Xi — > vqq. The relative branching ratio of the two modes is a function of the parameters of 
the model. The mode with a charged lepton is the better mode for reconstructing the Xi 
mass, because there is no neutrino which escapes detection. We have verified || that the 
charged lepton mode has a significant branching fraction over all of the mSUGRA parameters 
space allowed by present experimental searches. 

3 Xi reconstruction in mSUGRA models 

3.1 Introduction 

The possibility to extract a Xi mass peak is determined by: 

• the number of produced Xi f° r which the decay Xi ~* Z-jet-j et with I = e, /i can be 
completely reconstructed in the detector; 

• the combinatorial background, both from SUSY events and from Standard Model 
processes. 

These factors depend on the SUSY model considered, in particular the xi mass, the squark 
and gluino masses, and the jet and lepton multiplicity in SUSY events. Our aim is to 
determine if the reconstruction can be performed for the full parameter space. For this reason 
we chose to perform a detailed study for a few representative points in parameter space, 
trying to span a range as broad as possible of the parameters listed above. A convenient 
choice are the models already studied in detail for R-conserving mSUGRA , which present 
very varied phenomenologies, and also give the possibility to cross-check the results of our 
analysis with previous detailed studies. We have therefore selected points 1, 4 and 5, and 
we have added a low mass point, which we call 7, in a region with low branching fraction 
(~ 27%) for the decay xi ~~ " ^qq' ■ The mSUGRA parameters of the studied points are 
given in Table [I]. 
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Point 


m 


mi/2 


A 


tan p 


sgn(/x) 


1 


400 


400 





2.0 


+ 


4 


800 


200 





10.0 


+ 


5 


200 


300 


300 


2.1 


+ 


7 


200 


200 





2.0 


+ 



Table 1: Parameters of the example mSUGRA points studied. 



3.2 Selection criteria for jets 

In order to extract a significant x? peak the selection criteria need to be optimized separately 
for each addressed model. At this level, the important parameter is not the purity of the 
obtained Xi peak, but its significance, and the possibility of precisely measuring the peak 
position. 

An initial event selection was first made in order to reject the Standard model background. 
The events were required to have: 

• at least 6 jets with p t > 15 GeV 

• at least one lepton with p t > 20 GeV in point 4 and point 7, and at least two leptons 
with p t > 10 GeV in point 1 and point 5. 

The reason for requiring two leptons in point 1 and 5 is that a stronger ti background rejection 
is needed, since the SUSY cross-section for these points is lower. These are preliminary cuts 
and will be further optimized in the final analysis ||. For the Xi reconstruction the jets with 
the highest momentum were not considered, since there are likely to come from the start 
of the decay chain. The 3rd to the 8th jet in decreasing momentum scale were considered 
in point 1 and 5 and the 5th to the 8th jets in point 4 and 7. The Xi was reconstructed 
starting from an identified electron, and the invariant mass was calculated for the electron- 
jetl-jet2 combinations such that Ai?(electron-jl)< 2, Ai?(electron-j2)< 2, Ai?(jl-j2)< 2 
with AR = a/ A?] 2 + A0 2 , where t] is the pseudorapidity, and <p the angle in the plane 
perpendicular to the beam axis. In figure [I] the invariant mass distribution obtained in this 
way is shown in solid lines for the four points. A statistically significant peak, corresponding 
to the X? mass is seen in all four plots. The hatched histograms superimposed on the plots 
are the ti background. This important background can easily be evaluated and subtracted by 
performing the Xi reconstruction procedure using the wrong flavor of the lepton as described 
in the 'wrong flavor' technique below. In the next section we concentrate on the more difficult 
combinatorics from SUSY events. 

3.3 Subtraction of combinatorial background 

In order to estimate the shape of the combinatorial background two different techniques were 
combined. 
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Figure 1: Invariant mass of the accepted rnijj combinations for the sum of the SUSY signal 
and the SMti background for the four mSUGRA points studied (solid line). The background 
from ti production is shown as a hatched histogram. The dashed lines show the evaluation 
of the combinatorial background performed with the 'combined' method described in the text. 
The statistics shown correspond to an integrated luminosity of 30 fb^ 1 . 
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In the 'wrong flavor' technique, the combinatorial background is evaluated using all of 
the lepton-jet-jet combinations satisfying the cuts for the Xi candidate, but using the leptons 
with the wrong flavor, i.e. electrons if the Xi decays into muons and vice versa. 

In the 'mixed events' technique, the idea is to construct lepton-jet-jet combinations and 
make sure that not all objects in such a combination come from the same Xi> i- e - the 
combination does not contribute to the signal. In order to make sure that this is the case 
we use two different events (respectively 1 and 2). The directions (r),4>) of the three objects 
(lepton, jet, jet) and the momenta of two of the objects (lepton-jet or jet-jet) are taken from 
one of the candidate Xi combinations in event 1. The momentum of the third object is taken 
from a combination in event 2. 

These two approaches have complementary strengths and weaknesses. The 'wrong flavor' 
approach correctly subtracts all of the combinatorial background not involving one of the 
actual leptons from x\ decay, whereas the 'mixed events' background approximately accounts 
for events where only two of the three particles come from a Xi decay. 

In order to combine the strengths of the two methods the background is estimated through 
a linear combination of the backgrounds calculated with the two previous methods. The rel- 
ative weight of the two components is a function both of the characteristics of the SUSY 
model and of the applied cuts. In each case it is therefore necessary to optimize the param- 
eters of the linear combination. A prescription which was found to give reasonable results 
is: 

• Normalize the backgrounds estimated with the wrong lepton technique and with mixing 
events technique to the measured mass spectrum separately. This normalization is done 
in a region outside of the Xi peak. 

• Evaluate the relative weight of the two different background components by performing 
a least-square fit to the observed mijj mass distribution outside of the peak. 

The relative weight factors obtained exhibit some sensitivity to the choice of the mass 
windows over which the least square fit is performed. The 'real' background shape is best re- 
produced when the lower limit of the mass window is as near as possible to the peak position. 
Thus the estimate of the peak purity is affected by a systematic uncertainty, whereas the 
peak position has been checked to be independent of the choice of the subtraction technique. 
The sum of the background calculated with this combined method and the ti background is 
shown for all the four points as dashed lines in figure |I[ In all four cases the distribution for 
the combinatorial background is nicely reproduced. 

We show in Figure |2| the results for Point 5 in order to illustrate the power of this method. 
In the upper plot the invariant mass distribution for all the accepted electron-jet-jet combi- 
nations in SUSY events (solid line) is shown. The combinatorial background evaluated with 
the 'combined' method is superimposed as a hatched histogram. The lower left plot shows 
the signal peak after background subtraction (full line); the actual signal from Xi decay is 
superimposed (dashed line). Likewise, the full line in the lower right plot is the background 
distribution estimate obtained with the 'combined method', and the dashed line is the mass 
distribution for actual background combinatorics. Distributions showing the same excellent 
agreement are obtained for all four considered models. 
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Figure 2: Point 5: Invariant mass of the mijj combinations for SUSY events. In the upper 
plot the mass distribution for the combinatorial background evaluated with the 'combined' 
method (hatched) is shown superimposed on the observed signal+ background distribution 
(solid line). In the lower plot the mass distributions for signal (left) and background (right) 
are shown separately. The actual signal/background distributions are shown as dashed lines. 
The solid lines are the corresponding distributions evaluated with the 'combined ' method. The 
statistics shown corresponds to an integrated luminosity of 30 fb~ 1 . 
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4 Conclusions 



We have performed a reconstruction of Xi f° r the decay Xi ~~ * for different scenarios 
within the mSUGRA model with an R— parity violating A' term. We have shown that it is 
possible to define simple cuts on the jet and lepton multiplicity and topology which allow 
to observe a statistically significant peak in the lepton-jet-jet invariant mass distribution, 
at the position of the Xi mass. We have discussed a method to subtract the combinatorial 
background. Using this subtraction the Xi mass and the number of reconstructed Xi s can 
be precisely measured. This technique works effectively for all four points considered in 
this analysis. The mass of the reconstructed Xi can be used as a starting point for the 
reconstruction of particles further up in the decay chain. 
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Supersymmetry with R parity violation at the LHC: discovery 
potential from single top production 



P. CHIAPPETTA, A. DEANDREA, E. NAGY, S. NEGRONI, G. POLESELLO 

and J.M. VIREY 



Abstract 

We study single top production through imparity violating Yukawa type couplings 
at the LHC. We consider all 2 — ► 2 partonic processes at tree-level, including interfer- 
ence terms. The calculated 2^2 partonic cross sections are implemented in PYTHIA 
to generate complete particle final states. The generated events are processed through 
a fast particle level simulation of the ATLAS detector. We take into account all im- 
portant SM backgrounds and study the signal-to-background ratio as a function of the 
initial partonic states, the exchanged sparticle mass and width, and of the value of the 
Yukawa couplings. 

The feasibility of single top quark production via squark and slepton exchanges to probe 
several combinations of R parity violating couplings at hadron colliders has been studied 
According to those studies, the LHC is better at probing the B violating couplings A" whereas 
the Tevatron and the LHC have a similar sensitivity to A' couplings. We perform a complete 
and detailed study including all signal channels using a Monte Carlo generator based on 
Pythia 6.1 0, taking into account all the backgrounds and including the ATLAS detector 
response using ATLFAST 2.0 [§. 

The .R-parity violating parts of the Lagrangian that contribute to single top production 

are: 

% = KjA^hK - Kd^A + + h.c (i) 

The superscript c corresponds to charge conjugation. There are altogether 27 and 9 \'^ k and 
X'lj k Yukawa couplings, respectively. The most suppressed couplings are A' m , X' 133 , A' 112 , A 113 
(see H for detailed up to date reviews of the existing bounds). In order to fix the kinematical 
variables, the reaction we consider is 

Uiipx) + dj{p 2 ) -> t(p 3 ) + b{p 4 ) , (2) 

the pk being the 4-momenta of the particles and the indices i and j refer to the generations 
of the u and d-tjpe quarks. 

For valence- valence (VV) or sea-sea (SS) subprocesses, the scalar slepton exchange in the 
w-channel is taken into account but appears to be suppressed within our assumptions about 
the A' couplings and sfermion masses. The down type squark exchange in the s-channel 
squared amplitude is dominant. 

Let us now consider the subprocesses involving valence-sea (VS) quarks. Concerning R 
parity violating terms, slepton exchange in the s-channel and down type squark exchange in 
the M-channel contribute. The dominant terms are the squared amplitude due to e exchange, 
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and for initial quarks of the same generation (i = j), the interference between W and d. 
The result is sensitive to the interference term only if the product of A" couplings is large 
(around 1CT 1 ). For subprocesses involving quarks of different generations in the initial state 
the situation is more complex and all amplitudes have to be taken into account. 

We have carried out a feasibility study to detect single top production through i?-parity 
violation at the LHC by measuring the Ivbb final state using the following procedure. 

First, we have implemented the partonic 2^2 cross sections in the PYTHIA event 
generator. Providing PYTHIA with the flavors and momenta of the initial partons using a 
given parton distribution function (p.d.f.)Q, complete final states including initial and final 
state radiations and hadronization are generated. 

The events were processed through ATLFAST to simulate the response of the ATLAS 
detector. In particular the energy of electrons, photons and muons was smeared according 
to the resolution of the relevant detector element in the pseudorapidity range \r]\ < 2.5. 
Finally, a simple fixed cone algorithm (of radius R = 0.4) was used to reconstruct the parton 
jets. The minimum transverse energy of a jet was set at 15 GeV. According to the expected 
b-tagging performance of the ATLAS detector || for low luminosity at the LHC we have 
assumed a 60% b-tag efficiency for a factor 100 of rejection against light jets. 

The same procedure was applied to the SM background with the exception that we used 
besides PYTHIA also the ONETOP || event generator. 

The integrated luminosity for one year at low luminosity at the LHC is taken to be 10 
fb _1 . In Table [I] we display the total cross section values for different initial parton flavors 
in the case of exchanged squarks of mass of 600 GeV and of -R-parity conserving width 
= 0.5 GeV. We took for all A" = 10 _1 , which yields a natural width of the squark which is 
smaller than the experimental resolution. Table f| contains the same information for slepton 
exchange ( A' = 10 _1 , for a slepton of mass of 250 GeV and a width of r# = 0.5 GeV). Other 
processes are not quoted because the small value of the limits of their couplings prevents 
their detection. 



Initial partons 


cd 


cs 


ub 


cb 


Exchanged particle 


s 


d 


s 


d 


s 


Couplings 


A 212 A 332 


A 212 A 331 


^132^332 


A 231 A 331 


A 232 A 332 


Cross section in pb 


3.98 


1.45 


5.01 


0.659 



Table 1: Total cross-section in pb for squark exchange in the s-channel for a squark of mass 
of 600 GeV assuming T R = 0.5 GeV. 

In order to study the dependence of the signal on the mass and the width of the exchanged 
particle we have fixed the couplings to 10 _1 and have chosen three different masses for the 
exchanged squarks: 300, 600 and 900 GeV, respectively. For each mass value we have chosen 
two different T^: 0.5 and 20 GeV, respectively. For the first case dominates, whereas in 
the last one, when Y tot rs T r , the single top-production cross section decreases by a factor 

* We have used the CTEQ3L p.d.f. 
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~ 10. We have considered here the ub parton initial state, since this has the highest cross 
section value. 

In order to study the dependence on the parton initial state we have fixed the mass of 
the exchanged squark to 600 GeV and its width with Tr = 0.5 GeV and varied the initial 
state according to the first line of Table |T[ 

Finally, for the exchanged sleptons we have studied only one case, namely the ud initial 
state with a mass and width of the exchanged slepton of 250 GeV and 0.5 GeV, respectively. 
In each case we have generated about 10 5 signal events. 

The irreducible backgrounds are single top production through a virtual W (noted W*), 
or through W-gluon fusion. VT-gluon fusion is the dominant process (for a detailed study 
see fTj). A Wbb final state can be obtained either in direct production or through Wt or 
ti production. Finally, the reducible background consists of W+nj events where two of the 
jets are misidentified as 6-jets. 

For the tb final state first we reconstruct the top quark. The top quark can be recon- 
structed from the W and from one of the 6-quarks in the final state, requiring that their 
invariant mass satisfy 150 < < 200 GeV. The W can be in turn reconstructed from 

either of the two decay channels: W — > ud, W — > Iv. Here we have considered only the 
latter case which gives a better signature due to the presence of a high p t lepton and missing 
energy. The former case suffers from multi-jets event backgrounds. As we have only one 
neutrino, its longitudinal momentum can be reconstructed by using the W and top mass 
constraints. The procedure used is the following : 

- we select events with two b-jets of p t > 40 GeV, with one lepton of p t > 25 GeV, with 
pjmiss > 35 Q e y anc } w ith a jet multiplicity < 3, 

- we reconstruct the longitudinal component (p z ) of the neutrino by requiring Mi v = My/. 
This leads to an equation with twofold ambiguity on p z . 

- More than 80% of the events have at least one solution for p z . In case of two solutions, we 
calculate for each of the two b-jets and we keep the p z that minimizes \M top — Mi„b\. 

- we keep only events where 150 < M\ v \> < 200 GeV. 

Next, the reconstructed top quark is combined with the b quark not taking part in the top 
reconstruction. In order to reduce the ti background to a manageable level, we need to apply 
a strong jet veto on the third jet by requiring that its p t should be < 20 GeV. 

The invariant mass distribution of the tb final state after the cuts described above is 
shown in Fig. 1. Once an indication for a signal is found, we count the number of signal 
(N s ) and background (iV&) events in an interval corresponding to 2 standard deviations 
around the signal peak for an integrated luminosity of 30 fb _1 . Then we rescale the signal 
peak by a factor a such that N s /y/Nb = 5. By definition the scale-factor a determines the 
limit of sensitivity for the lowest value of the A" (A') coupling we can test with the LHC: 

Kjk ■ ALn < o.oi • yfc 

In Table |2| we show the limits obtained for the combinations of A'^A^ for different 
masses and widths of the exchanged s-quark. Also shown are the current limits assuming 
a mass for rhf = 100 GeV, the number of signal and background events, as well as the 
experimentally observable widths of the peak (T exp ). 
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Fig. 1 - Invariant mass distribution of Ivbb for the signal and backgrounds (dashed his- 
togram) after three years at LHC at low luminosity after having applied the cuts. 
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2.36xl0^ 3 


1.21xl0~ 2 


4.10xl0- 3 


1.51 xl0~ 2 


6.09xl0- 3 


2.09X10- 2 



Table 2: Limits for the values of the A'^A^ Yukawa couplings for an integrated luminosity 
of 30 fb _1 . For the other quantities see the text. The current limit is 6.25xl0 _1 . 
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In Table |3] we compile the sensitivity limit of the bilinear combination of the different 
Yukawa couplings one can obtain after 3 years of LHC run with low luminosity, if the 
exchanged squark has a mass of 600 GeV. 
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Couplings 


A 212 A 332 


A 212 A 331 


A 132 A 332 
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N s 


660 
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N b 
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T exp (GeV) 


38.5 


31.3 


37.5 


40.1 


Limits on A" x A" 


4.26xl0" 3 


7.08xl0" 3 


4.1xl0- 3 


l.llxlO -2 



Table 3: Limits on the Yukawa couplings for an exchanged squark of mass 600 GeV assuming 
r# = 0.5 GeV, for an integrated luminosity of 30 fb _1 . Current limit is 6.25xl0 _1 . 

For the exchanged sleptons (cf Table |j) we have calculated the sensitivity limit of the 
bilinear combination of the different Yukawa couplings only for the most favorable case, 
i.e. for the ud partonic initial state. We obtain 4.63xl0~ 3 for the limits on A' llfc A^3 3 (in 
comparison with the limit of 2.8xl0 -3 obtained by Oakes et al.). 
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A 13k A k33 


A 23k A k33 


Cross section in pb 


7.05 


4.45 


2.31 


1.07 


2.64 


0.525 



Table 4: Total cross-section in pb for slepton exchange in the s-channel for a slepton of mass 
of 250 GeV assuming T R = 0.5 GeV. 

For those cases where the exchanged squark (slepton) can be discovered at the LHC we 
have made an estimate on the precision with which one can determine its mass. For this 
purpose, we have subtracted the background under the mass peak and fitted a Gaussian 
curve on the remaining signal. For the assumed value of the coupling constant, the error on 
the mass determination is dominated by the 1% systematic uncertainty on the jet energy 
scale in ATLAS @. 
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Probing Jl p couplings through indirect effects on Drell-Yan 

production at the LHC 

D. CHOUDHURY and R.M. GODBOLE 
Abstract 

In this note we analyze the sensitivity to JR, P couplings at the LHC through the 
measurement of dilepton pairs, taking as example. We show that by exploiting 
the differences in the kinematic distributions of the R p contributions from that of the 
SM case, we can do significantly better than the bounds obtained from low energy 
experiments. Further the possible LHC bound on the couplings does not scale as the 
mass of the squark. Similar analysis can also be done to assess the reach of the LHC 
to effects due to leptoquarks/compositeness etc. 

Within the Standard Model, electroweak gauge invariance ensures both lepton number 
and baryon number conservation (at least in the perturbative context). However, this is 
not so within the Minimal Supersymmetric Standard Model (MSSM). In most standard 
treatments, such interactions are avoided by the imposition of a discrete symmetry called im- 
parity, which implies a conserved multiplicative quantum number, R = (—iy B + L + s ^ where 
B is baryon number, L is lepton number, and S is spin Q. All ordinary particles are 
i?-parity even, while all superpartners are i?-parity odd. If i?-parity is conserved, super- 
partners must be produced in pairs and the lightest superpartner, or LSP, is absolutely 
stable. Unfortunately, -R-parity invariance of the SUSY Lagrangian is an ad hoc assumption 
and not derived from any known fundamental principle. Hence, it is of interest to consider 
i?-parity-violating extensions of the MSSM, especially since the experimental signatures of 
low energy supersymmetry would then be radically different. Curiously, R p interactions 
can improve the agreement between theory and precision electroweak measurements, and 
also offer explanations for experimental anomalies such as the high-Q 2 excess reported at 
HERA §. 

The most general R p terms in the superpotential consistent with Lorentz invariance, 
gauge symmetry, and supersymmetry aref] 

Wft = KjkLiLjEl + \' ijk LiQjD c k + KikU'Dpi (1) 

where L, and Qi are the SU (2) doublet lepton and quark fields and Ef, Uf, Df are the 
singlet superfields. The U DD couplings violate baryon number while each of the other two 
sets violate lepton number. As our aim is to explore dilepton final states in pp collisions, we 
shall explicitly forbid the UDD interactions H as an economical way to avoid unacceptably 
rapid proton decay. 

The remarkable agreement between low-energy experimental data and SM expectations 
imply quite severe bounds on the strength of many R p operators ||. Thus, within the 

*We do not consider here the bilinear terms that mix lepton and Higgs superfields M. 



65 



context of colliders it is natural to consider the production of supersymmetric particles to 
be governed mainly by gauge interactions with R p being important only in the subsequent 
decay 0. However, such analyses are neither very sensitive to the actual size of the R p 
coupling nor are they particularly useful for the case of heavy superpartners (when the 
pair-production cross sections are smaller). A complementary approach is afforded by 
processes like Drell-Yan production where the exchange of a heavy squark is governed by 
the relevant R p coupling. In this note, we refine the analysis of Ref. [f| in the context of the 
LHC, thereby extending its reach. 

Expanding the superfield components in ([!]), we obtain the interaction Lagrangian that 
connects quarks to leptons: 

£-LQD = Kjk l^iL^-R^jL — ^iL<ikRUjL + dj L d kR U iL — UjL<ikR^iL + dl R l/ iL dj L — d C kR e ihUjh j +h.C. 

(2) 

With only one A^ fe being nonzero, the simplest processes (with an observable final state) to 
which Clqd would contribute are 

• Uj Uj -> e~e+ (4 R ) 

• d k d k -> e~ej (u jL ) 

• Ujdj -> efvi (4 R ) 

• QjQj -> m'j (#iL, e iL ) 

with the particle in parentheses being exchanged in the t-channel. If we allow more than 
one R p coupling to be non-zero, more exotic final states would be possible. However, as 
simultaneous existence of more than one such coupling is disfavored from the data on flavor- 



changing neutral current processes [10], we shall not consider this possibility. Furthermore, 
the limits on LQD couplings of muons with the first-generation quarks found in nucleon 
targets are weaker than those for electrons. Hence, in this note, we shall restrict ourselves to 
perhaps the most optimistic case, namely the study of a fi + fi~ final state in the presence of 
a non-zero A' 211 . The corresponding details for the other cases will be presented elsewhere. 

The signature we focus on is a without any missing transverse momentum. To 

ensure detectability, we demand that 

p T (/i ± ) > 50GeV and < 3 . (3) 

The lowest order SM process leading to such a final state is the Drell-Yan mechanism i.e. 
qq — > (7*, Z*) — > A non-zero A' 211 induces an additional t-channel diagram and the 

expressions for the total cross sections can be found in Ref. 0. QCD corrections to the 
Drell-Yan process have been calculated JTT]| to the next-to-leading order and are a function 



of the cm. energy (y^s) of the collider, the structure functions used and the subprocess 
scale M. The dependence on M is marginal though and one may approximate it by a 
scale-independent .fT-factor 1.12 (for the CTEQ4M |12| structure functions that we use). 



The same factor approximates the K-factor for the dimuon mass distributions also for the 
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range of dimuon masses and rapidities that we use. In the absence of a calculation of the 
higher-order corrections to the R p contribution, we assume that i^-factor for the full theory 
is the same as that within the SM. 

To maximize the sensitivity, we need to consider the differential distributions especially 
since the event topology of the R p "signal" is quite different from that of the SM "back- 
ground". A convenient set of independent kinematical variables is given by the dimuon 
invariant mass M, the rapidity of the pair ?7 pair and the difference of the individual 

rapidities Ar/. 
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Figure 1: The phase space distributions for the Drell-Yan process at the LHC. The kinemat- 
ical cuts are as in [|. For the rapidity distributions, an additional cut (M > 500 GeV) has 
been imposed. In each case, the lower curve corresponds to the SM and the upper curve to 
A' 211 = 0.5, m sq = 800 GeV. 

It is obvious that the contribution due to a heavy squark exchange would be non-negligible 
only for M > m sq . Thus the relative deviations are pronounced only for large values of M 
(see Fig.^J). It has been demonstrated [pl| that the parton luminosities can be measured 



with high accuracy at the LHC. Hence a change in the absolute value of the cross-section 
can indeed be used to look at the sensitivity of the signal to R p couplings. As a matter of 
fact such an analysis fl4| of the differential distribution in the dimuon mass of the dimuon 



pair cross-section has been used by the DO collaboration to put limits on the Quark-Electron 
Compositeness scale, using the Tevatron data. Recall that the analysis here does include the 
effect of the higher order corrections on the SM contribution to the pairs. 
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Note that Arj(= r] + — ?7_), the difference of the lepton rapidities is directly related to 
the scattering angle in the subprocess center-of-mass frame. Hence one expects to see the 
difference between an s-channel and a t-channel process in this distribution. This is borne 
out strikingly in the second of Fig. p]. The symmetry about At] = is dictated by the fact 
that the initial state is symmetric. And finally, one considers ^ pa i r which is a measure of the 
Lorentz boost of the center-of-mass frame, or in other words, the mismatch of the quark and 
antiquark energies. 

To quantify our comparison of the differential distributions in the two cases (SM vs. 
Jtp), we devise a \ 2 test. We divide the M-r] pau i T -Ar] hypervolume into equal sized bins and 
compare the number of signal (iVS M+ ^) and background (N^ M ) events in each bin for a given 
integrated luminosity. We then define a x 2 test of discrimination 
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where e is a measure of the systematic errors, arising mainly from the uncertainties in 
luminosity measurement and quark densities. To be specific, we use a uniform grid of 
(20GeV, 0.15, 0.15) and choose two representative values for e, namely 5% and 15%. 
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Figure 2: The 95% C.L. exclusion curves that may be obtained at the LHC for two different 
values of integrated luminosity. In each case, the lower curve corresponds to a systematic 
uncertainty (see eq. 0j e = 0.05 while the upper corresponds to e = 0.15. The shaded region 
indicates the region excluded by low energy measurements. 

In Fig. |2|, we present the exclusion plots in the (m sq , A' 211 ) plane that can be achieved at 
the LHC for a given luminosity. The interpretation is simple. The area above the respective 
curves can be ruled out at 95% C.L. As expected, for e = 0.15 we do slightly worse than 
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for the case of e = 0.05. And similarly, we do better with larger luminosity. What is most 
interesting, however, is that we do significantly better than the bounds obtained from 
low-energy experiments, viz. 

A' 211 < 0.059 (lOOGeV/mjJ . 

The LHC bound clearly does not scale as the mass of the squark. This is easily understood 
from a perusal of the first of Figs. p]. For smaller squark masses, the R p contribution is also 
peaked at smaller values of M. However, the SM contribution is also much bigger at such 
M. Consequently, the relative deviations are smaller and so are the contributions to the 
X 2 '■ For larger squark masses, the R p contributions concentrate at larger M and hence the 
relative deviations are larger too. This is also the reason why our analysis shows much more 
sensitivity to the R p couplings than the analysis of Ref. || . 

It is clear that a similar analysis can also be performed to study the sensitivity of LHC 
experiments to effects of Leptoquarks/compositeness etc. which would give rise to four- 
fermion currents with unusual chiral structure. 
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Neutrino masses, R-parity violating supersymmetry and collider 

signals 

A.K. DATTA, B. MUKHOPADHYAYA, S. ROY and F. VISSANI 

In a scenario with R-parity violating supersymmetry one can have masses for the neu- 
trinos. The seed of this lies in the superpotential, where, upon admission of lepton number 
violation, one can write the following terms over and above those of the minimal SUSY 
standard model (MSSM): 

W n = XijkhLjEl + K 3k UQPl + KjkU'Dpi + (1) 

If one considers the framework where R-parity is violated through a term of the type 
€iLiH2, then one neutrino acquires a mass at the tree level through mixing with neutralinos. 
We have argued that this may give rise to the mass splitting corresponding to the v il -^v T 
oscillation, which provides an explanation of the atmospheric muon neutrino deficit recorded 
by the superkamiokande (SK) experiment [Q. The other two neutrinos remain massless at 
the tree level. However, the A-and A -type terms can at the same time induce mass terms for 
all the three generations at the one-loop level, which can be responsible for a smaller mass 
splitting between the the first two neutrino generation, thereby explaining the solar neutrino 
problem. 

The relevant parameters in this scenario, in addition to those contained in the MSSM, 
are the e^'s mentioned above and vacuum expectation values (vev) for the scalar neutrinos, 
which are unavoidable consequences of the former. These can be lumped into one 'basis- 
independent' parameter, namely, the vev of the sneutrino corresponding to the state which 
acquires the tree-level mass. The SK results restrict this vev to be of the order of 100 keV. It 
is also possible to constrain the parameter space of the soft bilinear SUSY breaking terms of 
this theory from large angle mixing and neutrino mass-squared difference as demanded by the 
atmospheric data. In addition, non-trivial restrictions on the hierarchy of the e-parameters 
of different flavors or between the e's and /i follow if we demand that flavor-changing neutral 
currents be suppressed ||. Moreover, if one wants to solve the solar neutrino problem in this 
scenario, either through MSW solution or vacuum oscillation, then the A and A couplings 
also get restricted to values < 10~ 5 for the MSW solution, and to much smaller values for 
the 'just-so' solution. 

Apart from the neutrino-neutralino mixing mentioned above, such a scenario induces 
mixing in the char gino- charged lepton sector as well. These mixings result in some couplings 
like y^rW , x\vZ etc., which are characteristic of bilinear R-parity violation. 

If the lightest neutralino (xT) is heavier than the Z, then we have the following additional 
decay channels of Xi : 

xl^ViZ, l = e,n,T (2) 

Xi-W, l = e,fi,r (3) 

This kind of decay will not be possible with only the trilinear R-violating interactions. 
Also, if the A and A terms have to be responsible for loop-induced solutions to the solar 
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neutrino problem, then they are restricted to such values that three-body decays of the 
lightest neutralino triggered by them are dominated by the two-body decays as and when 
they are allowed. In such a case, demanding maximal mixing between the second and third 
generations (as required by the SK data), one would expect to get equal numbers of muons 
and taus in decays in the channel Xi ~~ * • This can be a remarkable feature in collider 
searches of the lightest neutralino. Another important observation is that, thanks to the 
small R-violating couplings necessitated by the small tree-level neutrino mass, the decay 
length of the lightest neutralino can be as large as 1 - 10 mm or so. This implies a decay 
gap that can be observed in collider experiments. 

Among other signals that have been investigated ||, the signal in the form of equal 
number of like sign fi and r events accompanied by like sign on-shell W's may be seen 
at the upgraded Tevatron, LHC or future generation e + e~ colliders when a pair of lightest 
neutralinos (xi) is produced directly or via cascades. Each neutralino may decay directly to 
yU or r (with BR ~ 40%) [[3J along with on-shell W's. We emphasize that this is a consequence 
of neutralino-neutrino mixing triggered exclusively by the bilinear R-parity violating terms 
in the superpotential. A very interesting possibility, essentially due to the Majorana nature 
of neutralinos, is that of having like-sign dilepton events accompanied by two W's of identical 
charge ||. Of course, for such a signal to be distinct, one requires the associated W's to 
be on-shell. In other words, such a signal is typical of a situation where Xi is heavier than 
the W. Otherwise, 3-body decays of the lightest neutralino take place, and these can have 
contributions from trilinear R-violating couplings as well. 

The interesting feature is that there is almost no Standard Model background for such 
a signal. Taming the MSSM background, too, is unlikely to pose any serious difficulty, so 
long as the W's can be identified |J. One may further exploit the large decay gap (~ few 
mm) P| in the lightest neutralino decay for strengthening the signal. 
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Introduction to GMSB phenomenology at TeV colliders 

S. AMBROSANIOQ 



Abstract 

A brief account is given of models with gauge-mediated supersymmetry breaking, as 
a general introduction for the following two contributions. 

1 Introduction to Gauge-Mediated SUSY Breaking 

Since no superpartners have been detected at collider experiments so far, supersymmetry 
(SUSY) cannot be an exact symmetry of Nature. The requirement of "soft" supersymmetry 
breaking |I| alone is not sufficient to reduce the free parameters to a number suitable for 
predictive phenomenological studies. Hence, motivated theoretical hypotheses on the na- 
ture of SUSY breaking and the mechanism through which it is transmitted to the visible 
sector of the theory [here assumed to be the one predicted by the minimal SUSY extension 
of the standard model (MSSM)] are highly desirable. If SUSY is broken at energies of the 
order of the Planck mass and the SUSY breaking sector communicates with the MSSM sec- 
tor through gravitational interactions only, one falls in the supergravity-inspired (SUGRA) 
scheme. The most recognized alternative to SUGRA is based instead on the hypothesis 
that SUSY breaking occurs at relatively low energy scales and is mediated mainly by gauge 
interactions (GMSB) A good theoretical reason to consider such a possibility is that 

it provides a natural, automatic suppression of the SUSY contributions to flavor-changing 
neutral currents and CP- violating processes. A pleasant consequence is that, at least in the 
simplest versions of GMSB, the MSSM spectrum and other observables depend on just a 
handful of parameters, typically 

M mess , N mcss , A, tan/3, sign(/x), (1) 

where M mess is the overall messenger scale; N mess is the so-called messenger index, parame- 
terizing the structure of the messenger sector; A is the universal soft SUSY breaking scale 
felt by the low-energy sector; tan j3 is the ratio of the vacuum expectation values of the two 
Higgs doublets; sign(/x) is the ambiguity left for the SUSY higgsino mass after conditions for 
correct electroweak symmetry breaking (EWSB) are imposed (see e.g. Refs. P-§|). 

The phenomenology of GMSB (and, more generally, of any theory with low-energy SUSY 
breaking) is characterized by the presence of a very light gravitino G ||, 

F ( VF V 

m 3/2 = = -= ~ — 2.37 eV, (2) 

1 ° V3M' P y 100 TeV J K J 
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where vF is the fundamental scale of SUSY breaking, 100 TeV is a typical value for it, 
and M' P = 2.44 x 10 18 GeV is the reduced Planck mass. Hence, the G is always the lightest 
SUSY particle (LSP) in these theories. If -R-parity is assumed to be conserved, any produced 
MSSM particle will finally decay into the gravitino. Depending on \/~F, the interactions of 
the gravitino, although much weaker than gauge and Yukawa interactions, can still be strong 
enough to be of relevance for collider physics. As a result, in most cases the last step of 
any SUSY decay chain is the decay of the next-to-lightest SUSY particle (NLSP), which can 
occur outside or inside a typical detector or even close to the interaction point. The pattern 
of the resulting spectacular signatures is determined by the identity of the NLSP and its 
lifetime before decaying into the G, 

1 ( VF \ f wnlsp r 5 
CTNLSP " 10OB (JOOTWJ 1 100 GeW ' [6) 

where B is a number of order unity depending on the nature of the NLSP. 

The identity of the NLSP [or, to be more precise, the identity of the sparticle(s) having 
a large branching ratio (BR) for decaying into the gravitino and the relevant SM partner] 
determines four main scenarios giving rise to qualitatively different phenomenology: 

Neutralino NLSP scenario: Occurs whenever < (m^ —m T ). Here typically a decay 
of the Ni to G7 is the final step of decay chains following any SUSY production 
process. As a consequence, the main inclusive signature at colliders is prompt or 
displaced photon pairs + X + missing energy. Ni decays to GZ° and other minor 
channels may also be relevant at TeV colliders. 

Stau NLSP scenario: Defined by m fl < Minfm^, mi ] — m T , features fi — > Gr decays, 
producing r pairs or charged semi-stable fi tracks or decay kinks + X + missing energy. 
Here and in the following, I stands for e or \x. 

Slepton co-NLSP scenario: When rug < Minfm^, + m T ], £r — > Gi decays are also 
open with large BR. In addition to the signatures of the stau NLSP scenario, one also 
gets £ + £~ pairs or £r tracks or decay kinks. 

Neutralino-stau co-NLSP scenario: If \m fl — < m T and < Tn^ R , both signa- 
tures of the neutralino NLSP and stau NLSP scenario are present at the same time, 
since N% <-> f\ 2-body decays are not allowed by phase space. 

Note that in the GMSB parameter space the relation mi > always holds. Also, 
one should keep in mind that the classification above is only valid as an indicative scheme 
in the limit m e , — > 0, neglecting also those cases where a fine-tuned choice of \[F and 
the sparticle masses may give rise to competition between phase-space suppressed decay 
channels from one ordinary sparticle to another and sparticle decays to the gravitino [pT0f| . 

In this report, two important aspects of the GMSB phenomenology at TeV colliders will 
be treated: 
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(A) The consequences of the GMSB hypothesis on the light Higgs spectrum using the most 
accurate tools available today for model generation and calculation; 



(B) Studies and possible measurements at the LHC with the ATLAS detector in the stau 
NLSP or slepton co-NLSP scenarios, with focus on determining the fundamental SUSY 
breaking scale y/F. 

For this purpose, we generated about 30000 GMSB models under well defined hypotheses, 



using the program SUSYFIRE 11 , as described in the following section. 



2 GMSB Models 

In the GMSB framework, the pattern of the MSSM spectrum is simple, as all sparticle masses 
are generated in the same way and scale approximately with a single parameter A, which 
sets the amount of soft SUSY breaking felt by the visible sector. As a consequence, scalar 
and gaugino masses are related to each other at a high energy scale, which is not the case 
in other SUSY frameworks, e.g. SUGRA. Also, it is possible to impose other conditions at a 
lower scale to achieve EWSB and further reduce the dimension of the parameter space. 

To build our GMSB models, we adopt the usual phenomenological approach, in particu- 
lar following Ref. ||, where problems relevant for GMSB physics at TeV colliders were also 
approached. We do not specify the origin of the SUSY higgsino mass /i, nor do we assume 
that the analogous soft SUSY breaking parameter Bfi vanishes at the messenger scale. In- 
stead, we impose correct EWSB to trade \x and B\i for M z and tan /?, leaving the sign of n 
undetermined. However, we are aware that to build a satisfactory GMSB model one should 
also solve the latter problem in a more fundamental way, perhaps by providing a dynamical 
mechanism to generate [i and Bfi, possibly with values of the same order of magnitude. This 
might be accomplished radiatively through some new interactions. However, in this case the 
other soft terms in the Higgs potential, namely m% i2 , will be also affected and this will in 
turn change the values of and B/j, coming from EWSB conditions [f§-|6|]. Within the study 
(A), we are currently considering some "non-minimal" possibilities for GMSB models that 
to some extent take this problem into account, and we are trying to assess the impact on 



the light Higgs mass. We do not treat this topic here, but refer to [12| for further details. 

To determine the MSSM spectrum and low-energy parameters, we solve the renormal- 
ization group (RG) evolution with the following boundary conditions at the M mess scale, 



M a = N mess Ag(-^-)^, (a = 1,2,3) 

VM mess y 47r 



m 2 



respectively for the gaugino and the scalar masses. In Eq. (|4|), g and / are the one- loop 
and two-loop functions whose exact expressions can be found e.g. in Ref. [/J, and C a are 
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the quadratic Casimir invariants for the scalar fields. As usual, the scalar trilinear couplings 
Af are assumed to vanish at the messenger scale, as suggested by the fact that they (and 
not their squares) are generated via gauge interactions with the messenger fields at the two 
loop-level only. 

To single out the interesting region of the GMSB parameter space, we proceed as follows. 
Barring the case where a neutralino is the NLSP and decays outside the detector (large \fF), 
the GMSB signatures are very spectacular and are generally free from SM backgrounds. 
Keeping this in mind and being interested in GMSB phenomenology at future TeV colliders, 
we consider only models where the NLSP mass is larger than 100 GeV, assuming that searches 
at LEP and the Tevatron, if unsuccessful, will in the end exclude a softer spectrum in most 
cases. We require that M mcss > 1.01A, to prevent an excess of fine-tuning of the messenger 
masses, and that the mass of the lightest messenger scalar be at least 10 TeV. We also impose 
M mess > Mqut exp(— 125/iV mess ), to ensure perturbativity of gauge interactions up to the 
GUT scale. Further, we do not consider models with M mess > 10 5 A. As a result of this and 
other constraints, the messenger index N mess , which we assume to be an integer independent 
of the gauge group, cannot be larger than 8. To prevent the top Yukawa coupling from 
blowing up below the GUT scale, we require tan/5 > 1.2 (and in some cases > 1.5). This 
is also motivated by the current bounds from SUSY Higgs searches at LEP II ||13|1 . Models 
with tan/3 > 55 (with a mild dependence on A) are forbidden by the EWSB requirement 
and typically fail to give m 2 A > 0. 

To calculate the NLSP lifetime relevant to our study (B), one needs to specify the value of 
the fundamental SUSY breaking scale VF on a model-by-model basis. Using perturbativity 
arguments, for each given set of GMSB parameters, it is possible to determine a lower bound 
according to Ref. 0, 



for collider physics, although some semi- qualitative cosmological arguments are sometimes 
evoked. 

In order to generate our model samples using SUSYFIRE, we used logarithmic steps for A 
(between about 45 TeV/iV mcss and about 220 TeV / y/ N mcss , which corresponds to excluding 
models with sparticle masses above ~ 4 TeV), M mcss /A (between 1.01 and 10 5 ) and tan/3 
(between 1.2 and about 60), subject to the constraints described above. SUSYFIRE starts 
from the values of SM particle masses and gauge couplings at the weak scale and then evolves 
up to the messenger scale through RGEs. At the messenger scale, it imposes the boundary 
conditions @ for the soft sparticle masses and then evolves the full set of RGEs back to the 
weak scale. The decoupling of each sparticle at the proper threshold is taken into account. 
Two-loop RGEs are used for gauge couplings, third generation Yukawa couplings and gaugino 
soft masses. The other RGEs are taken at the one-loop level. At the scale ^/m^m^, EWSB 
conditions are imposed by means of the one-loop effective potential approach, including 
corrections from stops, sbottoms and staus. The program then evolves up again to M mess 
and so on. Three or four iterations are usually enough to get a good approximation for the 
MSSM spectrum. 
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The light Higgs boson spectrum in GMSB models 

S. AMBROSANIO, S. HEINEMEYER and G. WEIGLEIN 



Abstract 

This study deals with the characteristics of the light Higgs boson spectrum allowed 
by the (minimal and non-minimal) GMSB framework. Today's most accurate GMSB 
model generation and two-loop Feynman-diagrammatic calculation of rrih have been 
combined. The Higgs masses are shown in dependence of various model parameters at 
the messenger and electroweak scales. In the minimal model, an upper limit on rrih of 
about 124 GeV is found for m t = 175 GeV. 



1 Introduction 

Within the MSSM, the masses of the CP-even neutral Higgs bosons are calculable in terms 
of the other low-energy parameters. The mass of the lightest Higgs boson, rrih, has been 
of particular interest, as it is bounded to be smaller than the Z° boson mass at the tree 
level. The one-loop results |]-§J for rrih have been supplemented in the last years with the 
leading two-loop corrections, performed in the renormalization group (RG) approach 
in the effective potential approach |7j and most recently in the Feynman-diagrammatic (FD) 
approach P|HTU|. The two-loop corrections have turned out to be sizeable. They can lower 
the one-loop results by up to 20%. These calculations predict an upper bound on rrih of 
about rrih < 130 GeV for an unconstrained MSSM with m t = 175 GeV and a common 
SUSY mass scale Msusy < 1 TeV. 



As discussed in the Introduction [pTJ| , the GMSB scenario provides a relatively simple 
set of constraints and thus constitutes a very predictive and readily testable realization of 
the MSSM. The main goal of the present analysis is to study the spectrum of the lightest 
neutral CP-even Higgs boson, rrih, within the GMSB framework. Particular emphasis is 
given to the maximal value of rrih achievable in GMSB after an exhaustive scanning of the 
parameter space. Our results are discussed in terms of the GMSB constraints on the low- 
energy parameters and compared to the cases of a SUGRA-inspired or an unconstrained 
MSSM. 



2 Calculation of rrih 

We employ the currently most accurate calculation to evaluate rrih, based on the FD approach 
as given in Refs. |8HlO|. The most important radiative corrections to rrih arise from the top 
and scalar top sector of the MSSM, with the input parameters m t , the masses of the scalar 
top quarks, m t ~ i; m t ~ 2 , and the t-mixing angle, 9+. Here we adopt the conventions introduced 
in Ref. ||. The complete diagrammatic one- loop result || has been combined with the 
dominant two-loop corrections of 0(aa s ) |||J and with the subdominant corrections of 
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m t ) GMSB models are generated with the program SUSYFIRE, according to the 



discussion of ||11|| . For this study, we consider only models with tan/3 > 1.5 (12| and rriA > 80 
GeV ||13|| . In addition, we always use m t = 175 GeV. A change of 1 GeV in m t translates 
roughly into a shift of 1 GeV (with the same sign) in rrih as well. Thus, changing m t affects 
our results on rrih in an easily predictable way. 

The results of the rrih calculation have been implemented in the Fortran code FeynHiggs 



14}| . This program has been combined with SUSYFIRE, which has been used to calculate the 



low energy parameters m 



ti- 



/i, M l7 M 2 , 



m 



a- 



for each of the ~30000 GMSB 



models generated. These have then been passed to FeynHiggs for the rrih evaluation in 
a coherent way. Indeed, we transfer the MS parameters in the SUSYFIRE output to on- 
shell parameters before feeding them into FeynHiggs. Compact expressions for the relevant 
transition formulas can be found in Refs. |15|,|T6 . 



Compared to an existing analysis in the GMSB framework [17]], we use a more complete 
evaluation of rrih- This leads in particular to smaller values of rrih for a given set of input 
parameters in our analysis. Also, in Ref. |17| although some GMSB scenarios with gener- 
alized messenger sectors were considered, the parameter space for the "minimal" case with 
a unique, integer messenger index N mcss — N± — N 2 — N 3 was not fully explored. Indeed, 
A was in most cases limited to values smaller than 100 TeV and M mcss was fixed to 10 5 
TeV. Furthermore, partly as a consequence of the above assumptions, the authors did not 
consider models with N mcss > 4, i.e. their requirements for perturbativity of the MSSM 
gauge couplings up to the GUT scale were stronger than ours. We will see in the following 
section that maximal rrih values in our analysis are instead obtained for larger values of the 
messenger scale and the messenger index. 



3 The Light Higgs Spectrum in GMSB 

In the following, we give some results in the form of scatter plots showing the pattern in 
GMSB for rrih, rriA as well as other low-energy parameters of relevance for the light Higgs 
spectrum. 

In Fig. |l](a), we show the dependence of rrih on tan/3, where only models with tan/3 > 
1.5, rriA > 80 GeV and ttinlsp > 100 GeV are considered, while m t is fixed to 175 GeV. 
The dependence is strong for small tan/3 < 10, while for larger tan/3 the increase of the 
lightest Higgs mass is rather mild. The maximum values for rrih — 124 GeV are achieved for 
tan j3 > 50. It should be noted that for very large tan j3 > 52, we also find a few models with 
relatively small rrih ^ 100 GeV. This is due to the fact that in this case EWSB conditions 



tend to drive rriA toward very small values [11]. This is made visible by the scatter plot in 



Fig. 0(b), where the pseudoscalar Higgs mass is shown as a function of tan/3. For such small 
values of rriA and for large tan/3, the relation rrih ~ tua holds. Thus small rrih values are 
quite natural in this region of the parameter space. On the other hand, one can see that 
extremely large values of rriA ^ 2 TeV can only be obtained for small or moderate tan /3 < 10 
GeV. A comparison between Fig |TJ(a) and (b) reveals that the largest m h values > 123 GeV 
correspond in GMSB to rriA values in the 300-800 GeV range. Indeed, it has been checked 
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Pseudoscalar Higgs Mass in GMSB 



(b) 




Figure 1: Scatter plots for the light scalar (a) and pseudoscalar (b) Higgs masses as functions 
of tan/3. Only GMSB models with tan/3 > 1.5, > 80 GeV and ttinlsp > 100 GeV are 
considered. 



that such large values are in general obtained in the FD calculation for 300 < vtla ^ 1000 
GeV, see Ref. f|. 

In Fig. ^j, we show the dependence of the lightest Higgs boson mass on the stop mixing 
parameter x top defined by 



NO 




Figure 2: The light CP-even Higgs boson mass is given as a function ofx top for tan (3 — 3, 50, 
rriA = 1000 GeV, a common soft SUSY breaking scale for the squarks, rrtq = 1000 GeV, and 
a gluino mass m~ g = 800 GeV. 



_ Aop — /Vtan/^ w h ere ms = J(m$ + m? )/2. (1) 
ms V 11 12 

For equal soft SUSY breaking parameters in the stop sector with the D-terms neglected, 
x top corresponds to the ratio X t /Ms of the off-diagonal and diagonal entries in the stop 
mixing matrix, see e.g. Ref. [jig] . 

Maximal values are obtained for x top ~ ±2, a minimum is reached around Xt op ~ 0. 
Thus, for large m-h values a large numerator in Eq. (fj) is required. From Fig. 0(a), one 
can see that in GMSB only negative values of A top are allowed at the electroweak scale, as 
a consequence of the fact that the trilinear couplings are negligible at the messenger scale. 
Due to the logarithmic dependence of m h on the stop masses, relatively large values of |-A top | 
are needed for large mn- In addition, large tan/3 is also required. From Fig. 0(b) one can 
check that this leads to values of x top ~ —0.95, which can only be achieved for positive \i. 
Fig. |](a) shows the dependence of A top on \i. Large values of |A top | are only reached for 
large \\x\ values. Therefore maximal h masses are obtained for relatively large and positive 
/i, as can be seen in Fig. f|(b)|] 

^In general, for large values of \/i\ and tan/3 the effects of the corrections from the b-b sector can become 
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Figure 3: Scatter plots of A top vs. ms, the mass scale appearing in the denominator of 
Eq. (a) and tan (3 vs. x top (b). 



All these arguments about the combination of low energy parameters needed for large 
m/j in GMSB are summarized in Tab. [1]. where we report the 10 models in our sample 
that give rise to the highest rrih values. Together with rrih, Tab. H shows the corresponding 
input GMSB parameters (Eq. ([I]) in [|TT|] ) as well as the values of the low energy parameters 
mentioned above. 

It is interesting to note that all the models shown in Tab. 1 feature a large messenger 
index and values of the messenger scale not far from the maximum we allowed while gener- 
ating GMSB models. We could not construct a single model with rrih ^ 122.5 GeV having 
-^Vmess < 6 or M mess < 10 5 TeV, for m t = 175 GeV. It is hence worth mentioning here that 
our choice of imposing M mcss /A < 10 5 =^> M m css <2x 10 10 GeV does not correspond to any 
solid theoretical prejudice. On the other hand it is true that M mcss > 3 x 10 s GeV always 
corresponds to gravitino masses larger than ~ 1 keV, due to Eqs. @ and ([5]) in ||11|| . The 
latter circumstance might be disfavored by cosmological arguments [0. A curious conse- 
quence is that the GMSB models with the highest rrih belong always to the stau NLSP or 
slepton co-NLSP scenarios. 

Note also that restricting ourselves to GMSB models with A < 100 TeV, M mess < 10 5 
TeV and N mcss < 4, we find a maximal rrih value of 122.2 GeV, for m t = 175 GeV and 



tan/? ~ 52. This is to be compared with the one-loop result of Ref. |T7|], m/,(max) = 131.7, 
for tan/3 around 30 (the assumed value of m t is not quoted). 

Values for rrih slightly larger than those we found here may also arise from non-minimal 



important, leading to a decrease in rrih. For the GMSB models under consideration, however, this is not the 
case as a consequence of the relatively large b masses. 
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contributions to the Higgs potential, in connection with a dynamical generation of \i and 
Bfi flll |. A treatment of this problem can be found in Ref. [ |19|1 . 



One should also keep in mind that our analysis still suffers from uncertainties due to 
unknown higher order corrections both in the RGEs for GMSB model generation and in the 
evaluation of mh from low energy parameters. A rough estimate of these effects leads to 
shifts in not larger than 3 to 5 GeV. 



4 Conclusions 

We conclude that in the minimal GMSB framework described above, values of rrih > 124.2 
GeV are not allowed for m t = 175 GeV. This is almost 6 GeV smaller than the maximum 
value for one can achieve in the MSSM without any constraints or assumptions about the 
structure of the theory at high energy scales P,PC|,|2~I]. On the other hand, the alternative 



mSUGRA framework allows values of that are ~ 3 GeV larger than in GMSB J[3| . This 
makes the GMSB scenario slightly easier to explore via Higgs boson search. This result 
was expected in the light of the rather strong GMSB requirements, such as the presence 
of a unique soft SUSY breaking scale, the relative heaviness of the squarks and the gluino 
compared to non-strongly interacting sparticles, and the fact that the soft SUSY breaking 
trilinear couplings Af get nonzero values at the electroweak scale only by RGE evolution. 
Nevertheless, once the whole parameter space is explored, it is not true that mGMSB gives 
rise to nih values that are considerably smaller than in mSUGRA. Even smaller differences 
in the maximal rrih might be present when considering non-minimal, complex messenger 
sectors [|17| or additional contributions to the Higgs potential [19|,[22|]. In any case, as for 



mSUGRA, current LEP II or Tevatron data on Higgs boson searches are far from excluding 
mGMSB, and the upgraded Tevatron and the LHC will certainly be needed to deeply test 
any realistic SUSY model. 
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A 
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tan j3 


m A 
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GeV 




10 6 TeV 


TeV 


GeV 




GeV 


GeV 


GeV 


GeV 




Al 
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7 


1.00 
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-0.90 


A2 
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8 
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Table 1: The 10 GMSB models giving rise to the highest values in our sample. For 
each model, together with the light Higgs mass, we show the values of the GMSB input 
parameters and other low energy parameters of interest for calculating . 
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Measuring the SUSY breaking scale at the LHC in the slepton 
NLSP scenario of GMSB models 

S. AMBROSANIO, B. MELE, S. PETRARCA, G. POLESELLO and 

A. RIMOLDI 

Abstract 

This study is focused on the measurement of the fundamental SUSY breaking scale 
y/F at the LHC in the GMSB scenario where a stau is the next-to-lightest SUSY par- 
ticle (NLSP) and decays into a gravitino with ctnlsp in the range 0.5 m to 1 km. This 
implies the measurement of mass and lifetime of long lived sleptons. The identifica- 
tion is performed by determining the time of flight in the ATLAS muon chambers. 
Accessible range and precision on y/F achievable with a counting method are assessed. 

1 Introduction 

The fundamental scale of SUSY breaking y/F is perhaps the most important quantity to 
determine from phenomenology in a SUSY theory. In the mSUGRA framework, the gravitino 
mass sets the scale of the soft SUSY breaking masses in the MSSM (~ 0.1 — 1 TeV), so that 
y/F is typically large ~ 10 10 ~ n GeV (Eq. (§) in 0). As a consequence, the interactions 
of the G with the other MSSM particles ~ F~ l are too weak for the gravitino to be of 
relevance in collider physics and there is no direct way to access y/F experimentally. In 
GMSB theories, the situation is completely different. The soft SUSY breaking scale of the 
MSSM and the sparticle masses are set by gauge interactions between the messenger and 
low energy sectors to be ~ asM^ (Eq. (§j) in 0), so that typical A values are ~ 10 — 100 
TeV. On the other hand, y/F is subject to the lower bound (|5|) in |IJ only, which tells us that 
values well below 10 10 GeV and even as low as several tens of TeV are perfectly reasonable. 
The G is in this case the LSP and its interactions are strong enough to allow NLSP decays to 
the G inside a typical detector size. The latter circumstance gives us a chance for extracting 
y/F experimentally through a measurement of the NLSP mass and lifetime (Eq. (|3]) in ]![). 

Furthermore, the possibility of determining y/F with good precision opens a window 
on the physics of the SUSY breaking ("secluded") sector and the way this SUSY breaking 
is transmitted to the messenger sector. Indeed, the characteristic scale of SUSY breaking 
felt by the messengers (and hence the MSSM sector) given by \J F mess in Eq. (H) can also 
be determined once the MSSM spectrum is known. By comparing the measured values 
of y/F and yj F mess it might well be possible to get information on the way the secluded 
and messenger sectors communicate with each other. For instance, if it turns out that 
y/F\^^ <C y/F, then it is very likely that the communication occurs radiatively and the ratio 
\J F mcss / F is given by some loop factor. On the contrary, if the communication occurs via a 
direct interaction, this ratio is just given by a Yukawa-type coupling constant, with values 
& I, see Refs. @,§. 
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An experimental method to determine y~F at a TeV scale e + e~ collider through the 
measurement of the NLSP mass and lifetime was presented in Ref. [|J, in the neutralino 
NLSP scenario. Here, we are concerned with the same problem, but at a hadron collider, 
the LHC, and in the stau NLSP or slepton co-NLSP scenarios. These scenarios provide a 
great opportunity at the LHC, since the characteristic signatures with semi-stable charged 
tracks are muon-like, but come from massive sleptons with (3 significantly smaller than 1. 
In particular, we perform our simulations in the ATLAS muon detector, whose large size 
and excellent time resolution || allow a precision measurement of the slepton time of flight 
from the production vertex out to the muon chambers and hence of the slepton velocity. 
Moreover, in the stau NLSP or slepton co-NLSP scenarios, the knowledge of the NLSP mass 
and lifetime is sufficient to determine y/F, since the factor B in Eq.(^) of |IJ is exactly equal 
to 1. This is not the case in the neutralino NLSP scenario where B depends at least on 
the neutralino physical composition and more information and measurements are needed for 
extracting a precise value of \/F. 

2 Choice of the Sample Models and Event Simulation 

The two main parameters affecting the experimental measurement at the LHC of the slepton 
NLSP properties are the slepton mass and momentum distribution. Indeed, at a hadron 
collider most of the NLSPs come from squark and gluino production followed by cascade 
decays. Thus, the momentum distribution is in general a function of the whole MSSM 
spectrum. However, one can approximately assume that most of the information on the 
NLSP momentum distribution is provided by the squark mass scale mq only (in the stau 
NLSP scenario or slepton co-NLSP scenarios of GMSB, one generally finds rrig > m,-). To 
perform detailed simulations, we select a representative set of GMSB models generated by 
SUSYFIRE. We limit ourselves to models with m NL sp > 100 GeV, motivated by the discussion 
in [|1|, and rrig < 2 TeV, in order to yield an adequate event statistics after a three-year low- 
luminosity run (corresponding to 30 fb" 1 ) at the LHC. Within these ranges, we choose eight 
extreme points (four in the stau NLSP scenario and four in the slepton co-NLSP scenario) 
allowed by GMSB in the (h-nlsp; m q) plane, in order to cover the various possibilities. 

In Tab. |T], we list the input GMSB parameters we used, while in Tab. |2| we report the 
corresponding values of the stau mass, the squark mass scale and the gluino mass. The 
"NLSP" column indicates whether the model belongs to the stau NLSP or slepton co-NLSP 
scenario. The last column gives the total cross section in pb for producing any pairs of SUSY 
particles at the LHC. 

For each model, the events were generated with the ISAJET Monte Carlo || that in- 
corporates the calculation of the SUSY mass spectrum and branching fraction using the 
GMSB parameters as input. We have checked that for the eight model points considered the 
sparticle masses calculated with ISAJET are in good agreement with the output of SUSYFIRE. 

The generated events were then passed through ATLFAST |7J , a fast particle- level simula- 
tion of the ATLAS detector. The ATLFAST package, however, was only used to evaluate the 
efficiency of the calorimetric trigger that selects the GMSB events. The detailed response of 
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3 
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B3 
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5 
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+ 


B4 
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+ 


B6 


2.31xl0 5 


3 


65.2 


1.83 




B7 


7.57xl0 5 


3 


104 


8.54 
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Table 1: Input parameters of the sample GMSB models chosen for our study. 
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£ 
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B5 
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B6 
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£ 
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B7 


399.2 


f 
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B8 


302.9 


£ 
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2430 
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Table 2: Features of the sample GMSB model points studied. 



the detector to the slepton NLSP has been parameterized for this work using the results of 
a full simulation study, as described in the next section. 

3 Slepton detection 

The experimental signatures of heavy long-lived charged particles at a hadron collider have 
already been studied both in the framework of GMSB and in more general scenarios p|-|il"|1 . 
The two main observables one can use to separate these particles from muons are the high 
specific ionization and the time of flight in the detector. 

We concentrate here on the measurement of the time of flight, made possible by the 
timing precision ( < 1 ns) and the size of the ATLAS muon spectrometer. 

It was demonstrated with a full simulation of the ATLAS muon detector |I2] that the 
j3 of a particle can be measured with a resolution that can be approximately parameterized 
as <j(/3)//3 2 = 0.028. The resolution on the transverse momentum measurement for heavy 
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particles is found to be comparable to the one expected for muons. We have therefore 
simulated the detector response to NLSP sleptons by smearing the slepton momentum and 
(3 according to the parameterizations in Ref. fll2| . 

An important issue is the online selection of the SUSY events. We have not made any 
attempt to evaluate whether the heavy sleptons can be selected using the muon trigger. For 
the event selection, we rely on the calorimetric E^ 1SS trigger, consisting in the requirement 
of at least one hadronic jet with p? > 50 GeV, and a transverse momentum imbalance 
calculated only from the energy deposition in the calorimeter larger than 50 GeV. We checked 
that this trigger has an efficiency in excess of 80% for all the considered models. 

A detailed discussion of the experimental assumptions underlying the results presented 
here is given in Ref. |13| . 



4 Event Selection and Slepton Mass Measurement 

In order to select a clean sample of sleptons, we apply the following requirements: 

• at least one hadronic jet with > 50 GeV and a calorimetric E™ ss > 50 GeV (trigger 
requirement); 

• at least one candidate slepton satisfying the following cuts: 

— \rj\ <2.4 to ensure that the particle is in the acceptance of the muon trigger 
chamber, and therefore both coordinates can be measured; 

— Aneas < 0.91, where /3 mea s is the (3 of the particle measured with the time of flight 
in the precision chambers; 

— The Pt of the slepton candidate, after the energy loss in the calorimeters has 
been taken into account, must be larger than 10 GeV, to ensure that the particle 
traverse all of the muon stations. 

If we consider an integrated luminosity of 30 fb" 1 , a number of events ranging from a 
few hundred for the models with 2 TeV squark mass scale to a few hundred thousand for a 
500 GeV mass scale survive these cuts and can be used for measuring the NLSP properties. 

From the measurements of the slepton momentum and of particle j3, the mass can be 

determined using the standard relation m = p ^ ■ For each value of (3 and momentum, the 



measurement error is known and it is given by the parameterizations in Ref. |T2j . Therefore, 
the most straightforward way of measuring the mass is just to use the weighted average of 
all the masses calculated with the above formula. 

In order to perform this calculation, the particle momentum is needed, which implies 
measuring the r\ coordinate. In fact, with the precision chambers only one can only measure 
the momentum components transverse to the beam axis. 

The measurement of the second coordinate must be provided by the trigger chambers, for 
which only a limited time window around the beam crossing is read out, therefore restricting 
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the (3 range where this measurement is available. Hence, we have evaluated the achieved 
measurement precision for two different (3 intervals: 0.6 < (3 < 0.91 and 0.8 < (3 < 0.91 
for the eight sample points. We found a statistical error well below the 0.1% level for those 
model points having rriq < 1300 GeV. Even for the three models (B5, B7, B8) with lower 
statistics (rriq ~ 2 TeV), the error stays below the 0.4% level. 

Many more details, tables and figures about this part of our study can be found in 
Ref. M. 



5 Slepton Lifetime Measurement 

The measurement of the NLSP lifetime at a high energy e + e~ collider was studied in detail 
in Ref. | J for the neutralino NLSP case. Similar to that study, the measurement of the 
slepton NLSP lifetime we are interested in here can be performed by exploiting the fact that 
two NLSPs are produced in each event. One can therefore select N\ events where a slepton 
is detected through the time-of-flight measurement described above, count the number of 
times N 2 when a second slepton is observed and use this information to measure the lifetime. 
Although in principle very simple, in practice this method requires an excellent control on 
all possible sources of inefficiency for detecting the second slepton. 

We give here the basis of the method, without mentioning the experimental details. We 
provide an estimate of the achievable statistical error for the models considered and a pa- 
rameterization of the effect on the lifetime measurement of a generic systematic uncertainty 
on the slepton efficiency. In case the sparticle spectrum and BRs can be measured from the 
SUSY events, as e.g. shown in Ref. |14[], an accurate simulation of all the SUSY produc- 
tion processes can be performed, and the results from this section are representative of the 
measurement precision achievable in a real experiment. 

Another method based on the same principles, but assuming minimal knowledge of the 
SUSY spectrum, is described in Ref. [[OJ, where a detailed estimate of the achievable sys- 
tematic precision is given. 

We define N% starting from the event sample defined by the cuts discussed in Sec. [|, 
with the additional requirement that, for a given value of the slepton lifetime, at least one 
of the produced sleptons decays at at a distance from the interaction vertex > 10 m, and is 
therefore reconstructed in the muon system. For the events thus selected, we define N 2 as 
the subsample where a second particle with a transverse momentum > 10 GeV is identified 
in the muon system. The search for the second particle should be as inclusive as possible, 
in order to minimize the corrections to the ratio. In particular, the cut f3 mcSuS < 0.91 is not 
applied, but particles with a mass measured from (3 and momentum incompatible with the 
measured slepton mass are rejected. This leaves a background of high momentum muons in 
the sample that can be statistically subtracted using the momentum distribution of electrons. 
The ratio 

No 
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Figure 1: The ratio R = N2/N1 defined in the text as a function of the slepton lifetime cr. 
Only the curves corresponding to the model points Bl, B5, B6, B8 are shown. 

is a function of the slepton lifetime. Its dependence on the NLSP lifetime cr in meters 
in shown in Fig. H for four among our eight sample models. The curves for the model 
points not shown are either very similar to one of the curves we show or are mostly included 
between the external curves corresponding to points Bl and B8, thus providing no essential 
additional information. Note that the curve for model 6 starts from cr = 2.5 m and not from 
cr = 50 cm, as for the other models. This is due to the large value of M mess (cfr. Tab. [1]), 
determining a minimum NLSP lifetime allowed by theory which is macroscopic in this case 
(Eqs. © and (§) in §). 

The probability for a particle of mass m, momentum p and proper lifetime r to travel a 
distance L before decaying is given by the expression 

P(L) = e - mL / pCT . (2) 

N2 is therefore a function of the momentum distribution of the slepton, which is deter- 
mined by the details of the SUSY spectrum. One therefore needs to be able to simulate the 
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1 10 10 1 10 10 

ct (m) 

Figure 2: Fractional error on the measurement of the slepton lifetime cr, for model sample 
points Bl to B4. We assume an integrated luminosity of 30 fb _1 . The curves are shown 
for three different assumptions on the fractional systematic error on the R measurement: 
statistical error only (full line), 1% systematic error (dashed line), 5% systematic error 
(dotted line). 

full SUSY cascade decays in order to construct the ct-R relationship. 
The statistical error on R can be evaluated as 

Relevant for the precision with which the SUSY breaking scale can be measured is instead 



92 



, Point 5 , Point 6 

1 m % i 1 1 nun i i nihil-- i 1 1 1 en i 1 1 mm i 1 i hulk i i m 



10 — 



-2 

10 b- 



b» 3 

o 10 
o 




I 1 1 



I L_LL 



-= 10 =- 



-= 10 =- 



10 10 

Point 7 



-1 



-2 



10 



#l I I I Mill 


ttt 








. " " y ~? 
• ' * '.S^ 

' * * Ks'^ 


HI I 


i i i 



10 



10 



10 



-= 10 =- 



-= 10 =- 




10 



±_LL 



I l_LL 



10 



10 



cx (m) 

Figure 3: The same as in Fig. but for the model sample points B5 to B8. 



the error on the measured cr. This can be extracted from the curves shown in Fig. |1] and 
can be evaluated as 



(Tier 



a(R)/ 



dR{cr) 
Oct 



(4) 



The measurement precision calculated according to this formula is shown in Figs. |2| and 
H] for the eight sample points, for an integrated luminosity of 30 fb -1 . The full line in the 
plots is the error on cr considering the statistical error on R only. The available statistics is 
a function of the strongly interacting sparticles' mass scale. Even if a precise R-ct relation 
can be built from the knowledge of the model details, there will be a systematic uncertainty 
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in the evaluation of the losses in N 2 , because of sleptons produced outside the rj acceptance, 
or absorbed in the calorimeters, or escaping the calorimeter with a transverse momentum 
below the cuts. The full study of these uncertainties is in progress. 
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Figure 4: Fractional error on the measurement of the SUSY breaking scale ■s/F for model 
sample points Bl to B4. We assume an integrated luminosity of 30 fb -1 . The curves are 
shown for the three different assumptions on the fractional systematic error used in Figs. 
and £| 

At this level, we just parameterize the systematic error as a term proportional to R, 
added in quadrature to the statistical error. We choose two values, 1%R and 5%R, and 
propagate the error to the cr measurement. The results are represented by the dashed and 
dotted lines in Figs. |] and ^ 
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For the models with squark mass scales up to 1200 GeV, assuming a 1% systematic error 
on the measured ratio, a precision better than 10% on the cr measurement can be obtained 
for lifetimes between 0.5-1 m and 50-80 m. If the systematic uncertainty grows to 5%, the 
10% precision can only be achieved in the range 1-10 m. If the mass scale goes up to 2 TeV, 
even considering a pure statistical error only, a 10% precision is not achievable. However a 
20% precision is possible over cr ranges between 5 and 100 m, assuming a 1% systematic 
error. 



Point 5 Point 6 




100 150 200 250 300 50 100 150 200 

VF/100 TeV 

Figure 5: The same as in Fig. but for model sample points B5 to B8. 

Note that the curves corresponding to the model points B2, B6 and B7 do not start from 
cr = 50 cm, but from the theoretical lower limit on cr of 1.8, 2.5 and 6.1 meters, respectively. 
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6 Determining the SUSY Breaking Scale y/F 



Using the measured values of cr and the NLSP mass, the SUSY breaking scale \J~F can 
be calculated from Eq.(^) in ||, where 3 = 1 for the case where the NLSP is a slepton. 
From simple error propagation, the fractional uncertainty on the \/F measurement can be 
obtained adding in quadrature one fourth of the fractional error in cr and five fourths of the 
fractional error on the slepton mass. 

In Figs. f| and [|, we show the fractional error on the y/~F measurement as a function 
of \f~F~ for our three different assumptions on the cr error. The uncertainty is dominated 
by cr for the higher part of the \J~F range and grows quickly when approaching the lower 
limit on y/~F. This is because very few sleptons survive and the statistical error on both 
mi and cr gets very large. If we assume a 1% systematic error on the ratio R from which 
cr is measured (dashed lines in Figs. |] and the error on \/F is better than 10% for 
1000 < VF < 4000 TeV for model points B1-B4 with higher statistics. For points B5-B8, 
in general one can explore a range of higher \[F values with a small relative error, essentially 
due to the heaviness of the decaying NLSP in these models. Note also that the theoretical 
lower limit (D in || on \fF is equal to about 1200, 1500, 3900, 8900 TeV respectively in 
model points B2, B5, B6, B7, while it stays well below 1000 TeV for the other models. 



7 Conclusions 

We have discussed a simple method to measure at the LHC with the ATLAS detector the 
fundamental SUSY breaking scale \f~F in the GMSB scenarios where a slepton is the NLSP 
and decays to the gravitino with a lifetime in the range 0.5 m < ctnlsp ~ 1 km. This method 
requires the measurement of the time of flight of long lived sleptons and is based on counting 
events with one or two identified NLSPs. It relies on the assumptions that a good knowledge 
of the MSSM sparticle spectrum and BRs can be extracted from the observation of the SUSY 
events and that the systematic error in evaluating the slepton losses can be kept below the 
few percent level. We performed detailed, particle level simulations for eight representative 
GMSB models, some of them being particularly hard due to low statistics. We found that a 
level of precision of a few 10's % on the SUSY breaking scale measurement can be achieved 
in significant parts of the 1000 < \[F < 30000 TeV range, for all models considered. More 
details as well as a full study of the systematics associated with this procedure and another 



less "model-dependent" method to measure VF is presented in detail in Ref. [13 
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Anomaly mediated SUSY breaking at the LHC 

F.E. PAIGE and J. WELLS 



Abstract 



Anomaly Mediated SUSY Breaking models are reviewed. Possible signatures at the 
LHC for one case of the minimal realistic model are examined. 

1 Introduction 

The signatures for SUSY at the LHC depend very much on the SUSY masses, which presum- 
ably result from spontaneous SUSY breaking. It is not possible to break SUSY spontaneously 
using just the MSSM fields; instead one must do so in a hidden sector and then communi- 
cate the breaking through some interaction. In supergravity models, the communication is 
through gravity. In gauge mediated models it is through gauge interactions; the gravitino is 
then very light and can play an important role. Simple examples of both have been discussed 
previously. A third possibility is that the hidden sector does not have the right structure 
to provide masses through either mechanism; then the leading contributions come from a 
combination of gravity and anomalies. This is known as Anomaly Mediated SUSY Breaking 
(AMSB), and it predicts a different pattern of masses and signatures. 

2 Anomaly-Mediated Supersymmetry Breaking 

In the supersymmetric standard model there exist AMSB contributions to the soft mass 
parameters that arise via the superconformal anomaly The effect can be understood 

by recognizing several important features of supersymmetric theories. First, supersymmetry 
breaking can be represented by a chiral superfield $ = 1 + m^/ 2 9 2 which also acts as a 
compensator for super- Weyl transformations. Treating $ as a spurion, one can transform a 
theory into a super-conformally invariant theory. Even if a theory is superconformal at the 
outset (i.e., no dimensionful couplings), the spurion $ is employed since the quantum field 
theory requires a regulator that implies scale dependence (Pauli-Villars mass, renormaliza- 
tion scale in dimensional reduction, etc.). To preserve scale invariance the renormalization 
scale parameter \i in a quantum theory then becomes /x/v$t$. It is the dependence of 
the regulator on $ that induces supersymmetry breaking contributions to the scalars and 
gauginos. 

The anomaly induced masses can be derived straightforwardly for the scalar masses. The 
Kahler kinetic terms depend on wave function renormalization as in the following superfield 
operator, 




(1) 



98 



Taylor expanding Z around /x and projecting out the FF* terms yields a supersymmetry 
breaking mass for the scalar field Q: 

2 1 d 2 In Z Q 2 1 ( 8jq 8jq \ 2 /„x 



Similar calculations can be done for the gauginos and the A terms: 

M % = - 9 4^m 3/2 = ^m 3/2l (3) 
2 dm fx Qi 

where the sum over a includes all fields associated with the Yukawa coupling y in the super- 
potential. 

There are several important characteristics of the AMSB spectrum to note. First, the 
equations for the supersymmetry breaking contributions are scale invariant. That is, the 
value of the soft masses at any scale is obtained by simply plugging in the gauge couplings 
and Yukawa couplings at that scale into the above formulas. Second, the masses are related 
to the gravitino mass by a one loop suppression. In AMSB Mj ~ m 3 / 2 ai/4:7r, whereas 
in SUGRA M, ~ ^3/2 • While the AMSB contributions are always present in a theory 
independent of how supersymmetry breaking is accomplished, they may be highly suppressed 
compared to standard hidden sector models. Therefore, for AMSB to be the primary source 
of scalar masses, one needs to assume or arrange that supersymmetry breaking is not directly 
communicated from a hidden sector. This can be accomplished, for example, by assuming 
supersymmetry breaking on a distant brane [[!]. Finally, the squared masses of the sleptons 
are negative (tachyonic) because (3 g > for U(l) and SU(2) gauge groups. This problem 
rules out the simplest AMSB model based solely on eqs. ||-[|. 

Given the tachyonic slepton problem, it might seem most rational to view AMSB as a 
good idea that did not quite work out. However, there are many reasons to reflect more 
carefully on AMSB. As already mentioned above, AMSB contributions to scalar masses are 
always present if supersymmetry is broken. Soft masses in the MSSM come for free, whereas 
in all other successful theories of supersymmetry breaking a communication mechanism must 
be detailed. In particle, hidden sector models require singlets to give the gauginos an accept- 
able mass. In AMSB, singlets are not necessary. Also, there may be small variations on the 
AMSB idea that can produce a realistic spectrum and can have important phenomenological 
consequences. This is our motivation for writing this note. 



3 Two realistic minimal models of AMSB: mAMSB 
and DAMSB 

As we discussed in the introduction, the pure AMSB model gives negative squared masses 
for the sleptons, thus breaking electromagnetic gauge invariance, so some additional contri- 
butions must be included. The simplest assumption that solves this problem is to add at 
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the GUT scale a single universal scalar mass m\ to all the sfermions' squared masses. We 
will call this model mAMSB. The description and many phenomenological implications of 
this model are given in Refs. @,|4|]. The parameters of the model after the usual radiative 
electroweak symmetry breaking are then 

m , WI3/2, tan/3, sgn/i = ±. 

This model has been implemented in ISAJET 7.48 |J; a pre-release version of ISA JET has 
been used to generate the events for this analysis. 

For this note the AMSB parameters were chosen to be 

m = 200 GeV, m 3 / 2 = 35 TeV, tan (3 = 3, sgn fi = + 

For this choice of parameters the slepton squared masses are positive at the weak scale, but 
they are still negative at the GUT scale. This means that charge and color might be broken 
(CCB) at high temperatures in the early universe. However, at these high energies there are 
also large finite temperature effects on the mass, which are positive (symmetry restoration 
occurs at higher T). In fact, a large class of SUSY models with CCB minima naturally fall 
into the correct SM minimum when you carefully follow the evolution of the theory from 
high T to today. If CCB minima are excluded at all scales, then the value of m must be 
substantially larger, so the sleptons must be quite heavy. 

The masses from ISAJET 7.48 for this point are listed in Table ||. The mass spectrum 
has some similarity to that for SUGRA Point 5 studied previously the gluino and 

squark masses are similar, and the decays x° ~~ * ^ an d X2 ~ * xV 1 are allowed. Thus, many 
of the techniques developed for Point 5 are applicable here. But there are also important 
differences. In particular, the xf is nearly degenerate with the Xi, n °t with the x®- The 
mass splitting between the lightest chargino and the lightest neutralino must be calculated as 
the difference between the lightest eigenvalues of the full one-loop neutralino and chargino 
mass matrices. The mass splitting is always above m n ±, thereby allowing the two-body 
decays Xi ~ ¥ Xx + 7r± HHH- Decay lifetimes of x ± are always less than 10 cm over mAMSB 
parameter space, and are often less than 1 cm. 

Another unique feature of the spectrum is the near degeneracy of the £l and Ir sleptons. 
The mass splitting is || 

mf L - m} R ~ 0.037 ^-m| cos 2/3 + M\ In 

There is no symmetry requiring this degeneracy, but rather it is an astonishing accident and 
prediction of the mAMSB model. 

It is instructive to compare the masses from ISAJET with those calculated in Ref. [[3] 
to provide weak-scale input to ISAJET. These masses are listed in the right hand side of 
Table |I[ Since the agreement is clearly adequate for the purposes of the present study, no 
attempt has been made to understand or resolve the differences. It is clear, however, that 
if SUSY is discovered at the LHC and if masses or combinations of masses are measured 
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Table 1: Masses of the SUSY particles, in GeV, for the mAMSB model with mo = 200 GeV, 
m 3/2 = 35TeV, tan/3 = 3, and sgn/i = + from ISAJET (left side) and from Ref. [Q] (right 
side) using the ISAJET sign conventions. 
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with the expected precision, then more work is needed to compare the LHC results with 
theoretical models in a sufficiently reliable way. 

Another variation on AMSB is deflected AMSB (DAMSB). The idea is based on Ref. 
who demonstrated that realistic sparticle spectrums with non-tachyonic sleptons can be 
induced if a light modulus field X (SM singlet) is coupled to heavy, non-singlet vector-like 
messenger fields \&j and ^ 

W mcss = X^X^^i. 

To ensure gauge coupling unification we identify ^ and as 5 + 5 representations of SU(5). 
When the messengers are integrated out at some scale M , the beta functions do not match 
the AMSB masses, and the masses are deflected from the AMSB renormalization group 
trajectory. The subsequent evolution of the masses below Mo induces positive mass squared 
for the sleptons, and a reasonable spectrum can result. Although there may be additional 
significant parameters associated with the generation of the \x and term in the model, 
we assume for this discussion that they do not affect the spectra of the MSSM fields. The 
values of \x and are then obtained by requiring that the conditions for EWSB work out 
properly. 

The parameters of DAMSB are 

m 3 / 2 , n, M Q , tan/3, sgn/x = ± 

where n is the number of 5 + 5 messenger multiplets, and M is the scale at which the mes- 
sengers are integrated out. Practically, the spectrum is obtained by imposing the boundary 
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Table 2: Masses of the SUSY particles, in GeV, for the DAM SB model with n = 5, M = 
10 15 GeV, and tan/? = 4 from Ref. 0. 
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conditions at M , and then using SUSY soft mass renormalization group equations to evolve 
these masses down to the weak scale. Expressions for the boundary conditions can be found 
in Refs. |5, 10], and details on how to generate the low-energy spectrum are given in [M] 



The resulting spectrum of superpartners is substantially different from that of mAMSB. The 
most characteristic feature of the DAMSB spectrum is the near proximity of all superpart- 
ner masses. In Table |] we show the spectrum of a model with n — 5, M = 10 15 GeV, and 



tan (3 = 4 as given in [ 10[ . The LSP is the lightest neutralino, which is a Higgsino. (Actually, 
the LSP is the fermionic component of the modulus X, but the decay of Xi to it is much 
greater than collider time scales.) All the gauginos and squarks are between 300 GeV and 
500 GeV, while the sleptons and higgsinos are a bit lighter (~ 150 GeV to ~ 250 GeV) in 
this case. 

In summary, we have outlined two interesting directions to pursue in modifying AMSB 
to make a realistic spectrum. The first direction we call mAMSB, and is constructed by 
adding a common scalar mass to the sfermions at the GUT scale to solve the negative 
squared slepton mass problem of pure AMSB. The other direction that we outlined is de- 
flected anomaly mediation that is based on throwing the scalar masses off the pure AMSB 
renormalization group trajectory by integrating out heavy messenger states coupled to a 
modulus. The spectra of the two approaches are significantly different, and we should ex- 
pect the LHC signatures to be different as well. In this note, we study the mAMSB carefully 
in a few observables to demonstrate how it is distinctive from other, standard approaches to 
supersymmetry breaking, such as mSUGRA and GMSB. 
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4 LHC studies of the example mAMSB model point 

We now turn to a study of the example mAMSB spectra presented in Table 1. A sample of 
10 5 signal events was generated; since the total signal cross section is 16 nb, this corresponds 
to an integrated LHC luminosity of 6 fb" 1 . All distributions shown in this note are normalized 
to 10 fb -1 , corresponding to one year at low luminosity at the LHC. Events were selected by 
requiring 

• At least four jets with p T > 100, 50, 50, 50 GeV; 

• E T > min(100GeV,0.2M cff ); 

• Transverse sphericity St > 0.2; 

• M cS > 600 GeV; 

where the "effective mass" M e s is given by the scalar sum of the missing Et and the p^'s of 
the four hardest jets, 

M cS = E T + p T>1 + p T:2 + p T ,3 + PT,4- 

Standard model backgrounds from gluon and light quark jets, ti, W + jets, Z + jets, and 
WW have also been generated, generally with much less equivalent luminosity. The M e g 
distributions for the signal and the sum of all backgrounds with all except the last cut are 
shown in Figure [|. The ISAJET IDENT codes for the SUSY events contributing to this plot 
are also shown. It is clear from this plot that the Standard Model backgrounds are small 
with these cuts, as would be expected from previous studies 

The mass distribution for £ + £~ pairs with the same and opposite flavor is shown in 
Figure The opposite-flavor distribution is small, and there is a clear endpoint in the 
same-flavor distribution at 




= 213.6, 215.3 GeV 



corresponding to the endpoints for the decays X2 — > ^lr^ ~^ Xi^ + ^~- This is similar to 
what is seen in SUGRA Point 5, but in that case only one slepton contributes. It is clear 
from the e + e~ + — e ± fi T dilepton distribution with finer bins shown in the same figure 
that the endpoints for £r and £l cannot be resolved with the expected ATLAS dilepton mass 
resolution. More work is needed to see if the presence of two different endpoints could be 
inferred from the shape of the edge of the dilepton distribution. 

Since the main source for X2 is Qr ~^ X2Q1 information on the squark masses can be 
obtained by combining the leptons from ~~ > ^ decays with one of the two hardest jets in 
the event, since the hardest jets are generally products of the squark decays. Figure ^| shows 
the distribution for the smaller of the two £ + £~j masses formed with the two leptons and 
each of the two hardest jets in the event. The dashed curve in this figure shows the same 
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Figure 1: Left: Effective mass distribution for signal (curve) and Standard Model background 
(shaded). Right: IS A JET IDENT codes for all produced particles contributing to M eff 
distribution after cuts. The dominant contributions are q~L (21-26), g (29), and q~R (41-46). 



distribution for M« > 175 GeV, for which the backgrounds are smaller. Both distributions 
should have endpoints at the kinematic limit for q~R — > x% ~^ ^ 
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652.9 GeV. 



Figure [| also shows the mass distribution formed with each of the two leptons combined 
with the jet that gives the smaller of the two ££j masses. This should have a 3-body endpoint 
at 

(Mf - M 2 o)(M 2 - M~ y 



M 2 )(M 2 
M% 



1/2 



605.4 GeV. 



The branching ratio for b\ — > is very small, so the same distributions with 6-tagged jets 
contain only a handful of events and cannot be used to determine the b\ mass. 

The decay chain q^ — > y^q — > £ LR £^q — > x\i + £~q also implies a lower limit on the ££q 
mass for a given limit on z = cos 9* or equivalently on the ££ mass. For z > (or equivalently 
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M|| (GeV) M|| (GeV) 

Figure 2: Mass distribution for opposite sign dileptons. Left: Distributions for same flavor 
signal (solid curve), opposite flavor signal (dashed curve), and Standard Model same flavor 
background (shaded). Right: e + e~ + — e ± /u, T distribution for signal on a finer scale. 



Ma > M$ ax /y/2) this lower limit is 

— MfM| + 3M 2 M 2 M 2 - M^Ml - M%M* - M^MfM 2 - 
M\MlMl + 3M 2 M 2 M 2 - M e 4 M 2 + (Mf — M 2 )x (6) 
y^M 4 + M 4 )(M 2 + M 2 ) 2 + 2M 1 2 M|(M| - 6M 2 M e 2 + M e 4 ) 
M^ n = 376.6 GeV 

where M g , M 2 , M e , and Mi are the (average) squark, (average) slepton, and xl masses. 
To determine this lower edge, the larger of the two It 'j masses formed from two opposite-sign 
leptons and one of the two hardest jets is plotted in Figure ^. An endpoint at about the 
right value can clearly be seen. 

The £ + £~, £ + £~q, £^q, and lower £ + £~q edges provide four constraints on the four masses 
involved. Since the cross sections are similar to those for SUGRA Point 5, we take the 
errors at high luminosity to be negligible on the £ + £~ edge, 1% on the £ + £~q and £^q upper 
edges, and 2% on the £ + £~q lower edge. Random masses were generated within ±50% of 
their nominal values, and the x 2 f° r the four measurements with these errors were used to 
determine the probability for each set of masses. The resulting distribution for the Xi mass, 
also shown in Figure has a width of ±11.7%, about the same as for Point 5; the errors for 
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Figure 3: Mass distribution for e + e~ + — e ^ events combined with one of the two 

hardest jets. Left: £ + £~j mass distribution (solid) and same with M u > 175 GeV (dashed). 
Right: £ ± j mass distribution for the one of the two hardest jets that gives the smaller ££j 
mass. 



the other masses are also comparable. Of course, the masses being measured in this case are 
different: for example the squark mass is the average of the q~n rather than the qi masses. 

The leptons from Xi ~ * Xi^ ±l/ are very soft. This implies that the rate for events with 
one or three leptons or for two leptons with opposite flavor are all suppressed. Figure |5] 
shows as a solid histogram the multiplicity of leptons with > 10 GeV and \r]\ < 2.5 for 
the AMSB signal with a veto on hadronic r decays. The same figure shows the distribution 
for a model with the same weak-scale mass parameters except that the gaugino masses Mi 
and M2 are interchanged. This model has a wino xt approximately degenerate with the X2 
rather than with the Xv Clearly the AMSB model has a much smaller rate for single leptons 
and a somewhat smaller rate for three leptons; these rates can be used to distinguish AMSB 
and SUGRA-like models. 

While the decay X2 — > Xih is kinematically allowed, the branching ratio is only about 
0.03%. Other sources of h in SUSY events are also quite small, so in contrast to SUGRA 
Point 5 there is no strong h — > bb signal. However, there is a fairly large branching ratio 
for g — > bb, ti with b — > Xi^, t xfb, giving two hard b jets and hence structure in the 
distribution. For this analysis b jets were tagged by assuming that any B hadron with 
Pt,b > 10 GeV and \t)b\ < 2 is tagged with an efficiency e# = 60%; the jet with the smallest 

R = V(A V ) 2 + (Acj)) 2 
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Figure 4: Left: Lower edge from larger ££j mass combining e + e + fi + fi — e ± /i =F with one 
of the two hardest jets. Right: Model-independent fit for xi mass. 



was then taken to be b jets. The two hardest jets generally come from the squarks. To 
reconstruct g — > bb one of the two hardest jets, tagged as a b, was combined with any 
remaining jet, also tagged as a b. In addition to the standard multijet and $t cuts, a 
cut M cS > 1200 GeV was made to reduce the Standard Model background. The resulting 
distributions for the b jet multiplicity and for the smallest bb dijet mass are shown in Figure || 
The dijet mass should have an endpoint at the kinematic limit for g —>■ bib — ► Xxbb, 



\ 



(M| — M~)(M~ — M? 



M- 



418.7 GeV. 



While the figure is roughly consistent with this, the endpoint is not very sharp; more work 
is needed to assign an error and to understand the high mass tail. There should also be a 
bt endpoint resulting from g — > tt , t — > X\b, with M-+ ~ M^o and essentially invisible. Of 
course m t has to be kept in the formula. This would be an apparent strong flavor violation 
in gluino decays and so quite characteristic of these models. Reconstructing the top is more 
complicated, so this has not yet been studied. 

The splitting between the Xi an d Xi 1S very small in AMSB models. ISAJET gives a 
splitting of 0.189 GeV for this point and cr = 2.8 cm, with the dominant decay being the 
two-body mode Xi ~^ Xi 7r v ^ a a virtual W. Ref. || gives a somewhat smaller value of [i 
and so a smaller splitting. The lifetime is of course quite sensitive to the exact splitting. 
Since the pion or electron is soft and so difficult to reconstruct, it seems better to look 
for the tail of long-lived winos. The signature is an isolated stiff track in a fraction of the 
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Figure 5: Lepton multiplicity with a hadronic r veto. Solid: AMSB model. Dashed: Same 
but with M 1 ^ M 2 . 



events that ends in the tracking volume and produces no signal in the calorimeter or muon 
system. Figure |7] shows the radial track length distribution in units of cr for winos with 
\t]\ < 1 and the (generated) momentum distribution for those with Rt > lOcr. Note that 
the ATLAS detector has three layers of pixels with very low occupancy at radii of 4, 11, and 
14 cm and four double layers of silicon strips between 30 and 50 cm. It seems likely that the 
background for tracks that end after the pixel layers would be small. 

It is instructive to compare this signature to that for GMSB models with an NLSP 
slepton. Both models predict long-lived charged particles with (3 < 1. In the GMSB models, 
two NLSP sleptons occur in every SUSY event, and they decay into a hard e's, /i's, or r's plus 
nearly massless G's. In the AMSB models, only a fraction of the SUSY events contain long- 
lived charged tracks, and these decay into a soft pion or electron plus an invisible particle. 
A detailed tracking simulation should be done for both cases. 

Acknowledgements: This work was supported in part by the U.S. Department of Energy 
under Contract DE-AC02-98CH10886. We also acknowledge the support of the Les Houches 
Physics Center, where part of this work was done. 

References 

[1] L. Randall and R. Sundrum, hep-th/9810155. 



108 




N b M bb (GeV) 



Figure 6: Left: Multiplicity of b jets for signal (curve) and Standard Model background 
(shaded). Right: Smallest mass for pairs of b jets for signal (curve) and Standard Model 
background (shaded). Both plots have M e g- > 1200 GeV in addition to the standard cuts and 
include a b tagging efficiency e& = 60%. 

[2] G.F. Giudice, M.A. Luty, H. Murayama and R. Rattazzi, JHEP 12, 027 (1998), hep- 
ph/9810442. 

[3] T. Gherghetta, G.F. Giudice, and J.D. Wells, hep-ph/9904378. 
[4] J. L. Feng and T. Moroi, hep-ph/9907319. 

[5] H. Baer, F.E. Paige, S.D. Protopescu and X. Tata, hep-ph/9810440. 

[6] I. Hinchliffe, F.E. Paige, M.D. Shapiro, J. Soderqvist, and W. Yao, Phys. Rev. D55, 
5520 (1997). 

[7] ATLAS Collaboration, ATLAS Detector and Physics Performance Technical Design 
Report, CERN/LHCC/99-15. 

[8] J.L. Feng, T. Moroi, L. Randall, M. Strassler and S. Su, Phys. Rev. Lett. 83, 1731 
(1999), hep-ph/9904250. 

[9] A. Pomarol and R. Rattazzi, JHEP 05, 013 (1999), hep-ph/9903448. 

[10] R. Rattazzi, A. Strumia and J.D. Wells, hep-ph/9912390. 



109 




R T /cx p T (GeV) 



Figure 7: Left: Radial track length distributions for Xi hi the barrel region, \rj\ < 1. Right 
p distribution for Xi with Rt > lOcr. 
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